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h Einleitung. 


Eine lineare homogene Differentialgleichung r-ter Ordnung 
yim + played + + p,(2)Y =( 
kann bekanntlich auf verschiedene Weise als ein System von Differentialgleichungen 


erster Ordnung geschrieben werden: 
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y,(z) u.’ y,(z)p (8), k= 1, aa 
oder, wenn man n Lösungen (y,„(2),...y,(2)), =1,...n auf einmal betrachtet 
(1) Y'(z) = Y(xz)P(x), 


wobei Y(2) = (y„(2)), Y’(2) = (y,(2)), Piz) = (pula), , k=1,. 
Zum Studium dieser Gleichung haben Volterra und ae zwei Operationen 
eingeführt, nämlich die Derivation 


(2) D,Y(2) = Y(a)Y(a) 
und deren Umkehrung, die Produktintegration 


(3) fa + Plt)dt) = im I (1 + Pla,)(a,,ı —4)); 


wobei der Grenzübergang in gleicher Art wie bei der Riemannschen Integraldefinition 
ausgeführt wird. 

Der Zweck der vorliegenden Arbeit !) ist einmal, die Theorie dieser zwei Operationen 
etwas weiter zu führen, u. a. durch Hervorhebung der Verbindung zwischen dem Produkt- 
integral und einer Reihenentwicklung für eine Lösung von (1), die zuerst von Peano [1] 
angegeben wurde und die nach einer einfachen Transformation die formell allgemeinere 
Reihe, zu der Schlesinger [1], S. 41 durch sukzessive Approximationen gelangt ist, enthält. 

Dann aber soll zweitens die so entwickelte Theorie verwertet werden bei der Lösung 
verschiedener klassischer Probleme der Theorie der Differentialgleichungen: der Ein- 
teilung der Integrale von Poincareschen Differentialgleichungen; dem Darstellungs- 
problem von Helge v. Koch [1], das dahingehend verallgemeinert wird, daß man die 
Lösung einer Gleichung sucht, deren Koeffizient eindeutig und regulär analytisch in 
einem beliebig doppelt zusammenhängenden Bereich ist; und schließlich der Diskussion der 
unwesentlichen Singularität. 

Bei der Lösung des ersten Problems werden einige wichtige Sätze von Perron 
durch direkte Abschätzung gewisser Produktintegrale bewiesen. 

Die Lösung des zweiten Problems beruht auf einem Kunstgriff, der im Prinzip von 
Frobenius [1] herrührt: Die inhomogene Gleichung 

(4) Y'(z;e)= Y(z; e)P(a) + =” 
hat eine Lösung, die bis auf den Faktor x* eindeutig im Bereich ist, und diese Lösung ist 
eine meromorphe Funktion des Parameters po, dessen Pole unabhängig von x sind, so daß 
die Residuen Partikularlösungen der homogenen Gleichung sind. 

Eine Gleichung mit einer unwesentlichen Singularität in x = 0 ist dadurch charak- 
terisiert, daß sie eine Lösung von der Form xz4V(x) besitzt, wobei V(x) regulär in 0 ist. 
Horn [2] hat gezeigt, daß sie rational in eine kanonische Gleichung 


(5) Ya) = Y(a)(2 + Fia)) 


transformiert werden kann, wo F(z) regulär in 0 ist und umgekehrt, daß eine solche 
Gleichung eine unwesentliche Singularität in 0 besitzt. Was nun A und B anbetrifft, 
so werden ihre Eigenwerte gleich sein, aber ihre Normalformen brauchen — im Gegensatz 
zu dem, was Volterra [3], S. 25—29 behauptet — nicht äquivalent zu sein; wieweit sie es 
sind, ist abhängig von F(x). Es gilt dagegen folgender Satz: Eine Gleichung mit einer un- 
wesentlichen Singularität kann rational in eine normale Gleichung von obenstehender Form 





ı) Einen wesentlichen Teil hiervon bilden die ersten zwei Kapitel meiner in dänischer Sprache geschriebenen 
Habilitationsschrift: Om Matrixregning og dens Anvendelse paa Differens- og Differentialligninger, Kopenhagen 1930. 
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transformiert werden, wobei F(x) in 0 regulär ist und die Eigenwerte von B keine ganz- 
zahligen Differenzen =# 0 haben; eine solche Gleichung hat eine Lösung von der Form x4V (x) 


mit A= B und außerdem V(0) =1. 


Kapitel 1. 
Hilfssätze über Matrizen. 


$ 1. Relativmatrizen. 


Eine rechteckige Matrix mit m Zeilen und n Spalten nennen wir von der Ordnung 
(m, n). Ist diese Matrix in pg rechteckige Felder eingeteilt, die also jedes für sich wieder 


Matrizen sind, 


(1,1) A= (Ar), Au= (m), 
so sagen wir, die ursprüngliche Matrix sei als Aelativmatrix von der Relativordnung (p, q) 
geschrieben. 


Ist eine andere Matrix B in entsprechender Weise eingeteilt, aber so, daß die Produkte 
A;; B;ı existieren, also A;; ebenso viele Spalten wie B;; Zeilen besitzt, dann kann man AB 
bilden durch Multiplikation der Relativmatrizen, so als ob deren Elemente Zahlen wären !): 

(1, 2) AB= (2ZA,,B,). 

Wenn A einen positiven Rang r hat, dann enthält sie wenigstens eine invertierbare 
Teilmatrix von der beliebigen Ordnung qg <r. Gewöhnlich kann man annehmen, daß 
diese Teilmatrix in der oberen linken Ecke steht, da man sie sonst durch Multiplikation 
mit zwei passend gewählten Vertauschungsmatrizen ?) dorthin bringen kann. Es ist 
zuweilen von Bedeutung, daß man, wenn A invertierbar ist, die invertierbare Teilmatrix 
innerhalb der ersten g Zeilen (Spalten) wählen kann; das erkennt man aus der Laplace- 
schen Entwicklung der Det. A nach den ersten g Zeilen (Spalten). In diesem Fall braucht 
A also nur mit einer Vertauschungsmatrix multipliziert zu werden. 

Die besagte Teilmatrix nennen wir A,, und teilen A ein in 


A„,ÄA 
1,3 ä | 11 a) 
Ay, Ay 
Danach kann man A in ein Produkt von drei relativen Elementarmatrizen zerlegen 
1 0\ /A, 0 1 A 
7 Tan) 
(4, 4) Ay 1 fi Ay} \0 1 
Dabei ist 4 
(1,5) e = An, Ay= An — Ayı An Ar, 
An = AnAn, Ay=AnAn. 


In dem speziellen Fall, daß A quadratisch und invertierbar ist, hat auch A,, diese 
Eigenschaften, und die Zerlegung von A kann auf andere Weise geschehen, nämlich 


4,48) (ll) 


wobei 
(1,5a) Ayı = Ay An. 
Es kann vorkommen, daß die Einteilung von A von der Form (A, Ass) 


oder (>) ist, aber man erkennt leicht, wie die Zerlegung in Faktoren dann modi- 


fiziert werden muß. 


1) Vgl. z.B. van der Waerden [1], S. 116. 
2) Eine Vertauschungsmatrix entsteht, wenn man in der Einheitsmatrix einige Zeilen oder Spalten unter- 


einander vertauscht (vgl. Kronecker [1], S. 369). 





gs 
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$ 2. Normalform. 
Zwei quadratische Matrizen A und B heißen äquivalent, wenn zwischen ihnen eine 
Relation von der Form | 
(2,1) A=ABAT" 
besteht, wobei A eine invertierbare Matrix ist. Nach Weierstraß [1] ist jede Matrix A 
einer Matrix P von besonders einfacher Struktur äquivalent. Diese Normalform kann 
als relativquadratische Matrix 


P, - - 0 
(2, 2) Pe. 
B: 
geschrieben werden, wo jedes P; von der Form 
00--:0 
(2, 3) mara[1 22100) (ev. = (0)) 
00..-.10 


ist. Die Zahlen o,, die man die Eigenwerte von A nennt, sind Nullstellen der charakteristi- 
schen Funktion Det. (A — oe) der gegebenen Matrix, und wenn P, von der Ordnung 4, ist, 
dann ist (o— o,)"' ein Elementarteiler dieser Determinante }). 

Da wir Matrizen von der Form © noch häufig brauchen werden, sei schon hier be- 
merkt: wenn ©, und ©, von der Ordnung m und n sind, während C = (c,,) von der Ordnung 
(m, n) ist, dann ist 








P Be En 208 0 
(2, 4) ecer— Be u 0 Ben 
ee ee Cı,n 0 0 
A „+1 Cn—u.n 0 0 
Besonders ist 
0 0- 0 0 
„fe o---000 Pe! 
(2,5) Oi=[1 0:-:0 0 0|,::,„ei” =|9 0...0p 9 =9 
SE Zn Eee wi 0 


$ 3. Permutable Matrizen. 


Über permutable Matrizen erwähnen wir hier zwei Sätze, von denen wir im Folgen- 
den noch Gebrauch machen werden. 

I. Die notwendige und hinreichende Bedingung dafür, daß alle Matrizen, die mit 
einer gegebenen Matrix permutabel sind, auch untereinander permutabel sind, ist, daß die 
Elementarteiler der gegebenen Matrix von der Form (o— 0,)",.... (o—0,)” und o,,... 0, 
verschieden sind. (Volterra [2], S. 11.) 

II. Die Eigenwerte zweier permutabler Matrizen A und B kann man so assozüieren 
O1» 915 ++ +3 0, 0, daß die Eigenwerte von «A + BB und AB, &o, + Bo,bzw. 0,0, i=1,...n 
sind. & und ß bezeichnen dabei beliebige Skalare. (Frobenius [3].) 


$ 4. Ungleichungen. Die Matrix E. 
Sind alle Elemente in A und Breell, so sollen durch die Relationen 
(4,1) AsBundA=zB 





1) Vgl. z.B. v. d. Waerden [1], $$109, 110. 
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die Ungleichungen 

(4, 2) ax S bu bzw. ar 2 bir 
ausgedrückt werden. 

Den absoluten Betrag einer Matrix definieren wir als die Matrix der absoluten Be- 
träge ihrer Elemente 


(4, 3) 14|= (las|). 

Es gelten dann augenscheinlich die folgenden Ungleichungen 

(4, 4) 1AI—I21s|A+B|s|]A|+[|2B]| 
und 

(4, 5) |AB|s|jA|-|B|. 


Das Gleichheitszeichen in der letzten Relation gilt nur in speziellen Fällen, z. B. wenn 
A und B positiv sind. 


Sind alle Elemente in A absolut genommen kleiner oder gleich —, so schreiben wir 
(4, 6) |AlSaE. 


Die Matrix E, deren Elemente alle gleich - sind und von der wir noch häufig Ge- 


brauch machen werden, besitzt mehrere einfache Eigenschaften. Zuerst ist 


(4, 7) E=E, »>4. 
Weiter ist 
1 %ı u er, 7 & uaR 0 
(4, 8) Bam al... )=E(.:::. ‚ =JLa; 
n o = 0 4 j=1 
und 
...0 m 
(4, 8a) AE=|------ E, &: =Sa,, 
0... ml 


woraus folgt 
(4, 9) EAE=saE, «=: Yan. 
i,k=1 
Bei Anwendung von Relativmatrizen kann es praktisch sein, E in folgender Weise 
zu zerlegen: 


En *  : Em 
(4, 10) N 
Emı a Emm 
wobei Ej1 - - - Emm quadratisch sind. Es gelten dann die Relationen 


(4,11) EyEr = ui .Eu, 
wobei w; die Ordnung von Z;; ist. 


$ 5. Matrixfolgen. 
Gibt es zu jedem positiven e unendlich viele », so daß 
(5, 1) |A,—Al<eE 
ist, dann sagen wir die Matrixfolge (A,) besitzt die Häufungsmatrix A. 
Die Folge heißt konvergent, wenn sie nur eine Häufungsmatrix besitzt, d.h. wenn 
die obige Ungleichung für jedes hinreichend große » erfüllt ist, oder, was dasselbe be- 
sagt, wenn jede Elementfolge konvergiert. 
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Das allgemeine Konvergenzprinzip gilt offenbar auch für Matrixfolgen: Besteht 
für jedes positive e und jedes hinreichend große # und » die Ungleichung 


(5, 2) |A„—A,| <eE, 
so ist die Folge (A,) konvergent. 


Weiterhin führen wir den Satz von Bolzano-Weierstrass an: Jede beschränkte - 


Matrixfolge besitzt mindestens eine Häufungsmatrix, und jede Häufungsmatrix läßt 
sich darstellen als Grenzwert einer konvergenten Teilfolge. 

Auf Grund der Ungleichungen (4,4) und (4,5) kann man leicht die folgenden 
Rechenregeln beweisen: 

(5, 3) lim (A, +3B,)=A+B, limA,B,=AB, 


wenn A,>A und B,>B. 
Existiert weiterhin lim A,', so hat man nach (5, 3) 


”r>n 


1 = lim A, - lim A,", 


d.h. A ist invertierbar und 
(5, 4) lim A,'= 4A". 


Umgekehrt: Ist A invertierbar, dann können wir 
A,=All—H,), |IN,|<mE 





setzen, wobei 7, > 0. Da die Reihe B5 H# die Majorantenreihe 1 —+ - nm E=1+ Ps - :E 


hat, ist sie konvergent, und zwar ist ihre Summe (1 — H,)"'. A, ist folglich invertierbar 
für jedes hinlänglich große » und 


(5, 5) Ar =M—H)"TAtT>4A". 


$ 6. Matrixfunktionen. 
Eine Matrix, deren Elemente Funktionen von einer oder mehreren Variablen sind, 
nennen wir eine Matrixfunktion. 


Indem man den Differentialquotienten und das Integral einer Matrix definiert 
als die Matrizen, deren Elemente Differentialquotienten, bzw. Integrale der Elemente 
der betrachteten Matrix sind, erkennt man, daß man mit Matrixfunktionen genau so 
wie mit gewöhnlichen Funktionen rechnen kann. Man muß dabei nur die oben ent- 
wickelten Rechenregeln beachten, und weiterhin bemerken, daß die Mittelwertsätze nur 
für die gewöhnlichen Funktionen gelten. 


Wir werden uns deshalb im Folgenden ohne weiteres der für gewöhnliche Funk- 
tıonen gebräuchlichen Terminologie: Stetigkeit, Differentiierbarkeit, totales Differential, 
analytisch, Cauchys Integral, analytische Fortsetzung usw. bedienen, um Eigenschaften 
der Matrixfunktionen auszudrücken. Es ist aber hierbei zu bemerken, daß einzelne Ele- 
mente Eigenschaften haben können, die der Matrixfunktion als solcher nicht zukommen. 
Sind z. B. die Elemente in F(x) = (f;,(2)) durch Potenzreihen definiert, so setzt man 
F(x) analytisch fort, durch simultane Fortsetzung der einzelnen Funktionen /,(x), und 
auf dieser Grundlage baut man die vollständige analytische Funktion F(x) auf. Es kann 
dabei geschehen, daß einzelne der Funktionen /,,(x) in einem größeren Bereich fortsetz- 
bar sind als F(x), ebenso wie eine Singularität für F(x) keineswegs eine Singularität 
für jedes einzelne Element in F(x) zu sein braucht. 
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$ 7. Analytische Funktionen von Matrizen. 


Ist f(x) eine analytische Funktion, so definieren wir nach einem Vorschlag von 
Frobenius [2] die analytische Funktion f(x) der Matrix X = (z;,.) und deren Differential- 
quotienten durch das Cauchysche Integral 

„ y! . 
(7,1) = Z; [OK — Kar. 


2ri 


Dabei ist über eine geschlossene Kurve zu integrieren, die die Eigenwcrte von X, aber 
keine singulären Punkte von /(x) umschließt. Man erkennt unmittelbar, daß das eine 
analytische Funktion der n? Veränderlichen (z,,, - - ., Zn) ist, und daß X ein regulärer 
Punkt ist, wenn keiner ihrer Eigenwerte mit einer Singularität von f(x) zusammenfällt. 
Aus der Definition von f(X) ergibt sich unmittelbar der Satz: 
I. Wenn A mit X permutabel ist, ist es auch mit f{X) permutabel. 


Ist Z== A”' XA die Normalform von X, so ist augenscheinlich 
(7, 2) H{X) = AZAT. 


Eine erste Vorstellung über die Struktur der Matrix /{Z) kann man sich leicht machen. 
Ist nämlich =, die Diagonalmatrix, die dieselbe Diagonale wie Z hat, und setzen wir 


(7,3) z=3,+®, 
so wird nach (2,5) eine gewisse Potenz von © gleich Null sein: 
e=0. 
Da =, und © permutabel sind, wird 


und folglich 
SE) 
(7,3) nz) = 3 Teer. 


(7, 2) ergänzt durch (7, 4) nennen wir die erste kanonische Form von f(X). Ihre besondere 
Bedeutung beruht darauf, daß f”(=,) eine Diagonalmatrix ist. Setzt man 


we. 
(7,5) n-( Pair ) 
ae... 


so wird nämlich 


LM)" 0 
C—3)" = | PDF TFTRRDETYTE 
0 a) 
und hieraus erkennt man unmittelbar, daß 
EP) ---0 
(7, 6) Po ee erde ) 
0 IE) 


Die zweite kanonische Form, die genaueren Aufschluß über die Struktur von f(Z) 
gibt, bekommt man, wenn = als eine relative Diagonalmatrix 
u ::: 


0 
=) +0 1&---0 
07,7) = -( Eee ) s=&+9=I1  ....... (ev. = (£;)) 
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geschrieben wird. Dann ist 


und hieraus folgt 





tz): :- 0 
(7, 8) ftz -( EEE ) 
Dr En 
Auf jedes Diagonalelement dieser Relativmatrix wenden wir nun (7,4) und (2,5) an: 
f&) 0 .+.+0 
a1 (E,) (si) f&) 0 
(7, 9) ftz:) = > as vi MT = A 
1] fe’ (e) 
= ii f&) 


und gewinnen damit einen Einblick in die Struktur der Matrix /(Z). 
Eine dritte kanonische Form!) erhält man durch Anwendung der Partialbruch- 
zerlegung von ((— X )=', deren Pole ja &),-- -, &m sind. Wird 





m Pi 
(7,10) ei a Ar 


in (7,4) mit » = 0 eingesetzt, so bekommt man offenbar 


m Pi 
f 


(7,14) 1x) = I IT x, 


i=1 ı=0 7 
Durch Vergleich mit der ersten kanonischen Form ergibt sich die interessante Formel 


RO EINE), = ME)OA". 
In gewissen Fällen kann die Normalform von f(Z) leicht angegeben werden. Ist 


nämlich die Ordnung e, der Matrix Z, positiv und f’(£,) +0, so entspricht f(Z,) nur 
ein Elementarteiler (o— f(£,))“, und es existiert daher eine Matrix B,, für die 


K=,) = B,(f(&,) + 9,)Bz" 


ist. Sind diese Bedingungen für alle Werte von : erfüllt, dann gibt es also eine Matrix 


(9:72) 
4 ; 
die durch 


(7, 13) ft=) = B(f(&,) + @)BT 
die Reduktion auf die Normalform leistet. Andernfalls bedarf die Bestimmung der 
Normalform einer eingehenderen Untersuchung, die wir aber hier nicht vornehmen 
wollen, da sie für das Folgende nicht nötig ist. 

Die Formel (7, 4) ist ein Spezialfall eines Satzes von H. B. Phillips [1]: 

Il. Sind A und X permutabel und o,,&,;:.-; 0, &, deren assoziierte Eigenwerte, so 
gilt die Entwicklung 


© +(») 
(7, 14) a) = SM) (x ar, 


wenn &, im Regularıtätskreis um o, für f(x) liegt. 
1) Bromwich [1]. 











ET 


N. 
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Als Integrationsweg in (7,1) mit »= 0 wählt man eine Kurve, die die Kreise 
|x—e,|=]|£,—.o,|, aber keine Singularität für f(x), umschließt. Die Normalform 
Fee 


der Matrix ((— A)"(X— A), deren Eigenwerte nach Frobenius BI sind (vgl. 


$3, II), ist absolut genommen 


Infolgedessen ist die Reihe 
E-K"=-L(C—A)UX— A) 


gleichmäßig konvergent, und damit ist der Satz bewiesen. 
Als Beispiel führen wir an, daß 


’) 


wenn f(x) regulär in O0 ist. Durch Wh Fortsetzung erkennt man, daß die hieraus 
folgenden Formeln 


Pem = (f(x) — /{0))E 
f{zE)E = f(z)E 


für jedes x gelten, für das /(x) regulär ist. 


(7, 15) 


$ 8. Exponentialfunktion und Logarithmus ?). 
Die Exponentialfunktion und der Logarithmus sind durch die Gleichung 
(8, 1) eu X 
untereinander verbunden, was man z.B. auf Grund der Potenzreihen 


. ee) zer 
e> == m 
7 


0 
g1—X)=— YT, Jal<t,i=h,...n 


beweisen kann. Für |X| <£Z, &£<1 hat man nämlich 


e-logi1-X) _ zeme (1 — X)) u 1 os ur xr. 
y! 


v„=( „=() nu=y 





Da die Koeffizienten a‘ alle positiv sind, kann leicht bewiesen werden, daß die Doppel- 
reihe absolut konvergent ist. Wir dürfen deshalb die Summationsordnung ändern, wo- 


durch wir 
ar) 


—lgg(1-X) — -37 5% 


„=(0 


bekommen, und hierin müssen die Koeffizienten von X” gleich 1 sein, was doch für ska- 
lares x der Fall ist. Folglich ist 

e-lsu-n - (1— X)", 
so daß (8, 1) erfüllt ist, wenn X in der Nähe von O liegt. Durch analytische Fortsetzung 
erkennt man, daß die Relation allgemeingültig ist. 





1) Diese Funktionen sind in gewisser Beziehung von J. v. Neumann [1] sorgfältig diskutiert worden. 
Journal für Mathematik. Bd. 171. Heft 2. 10 
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Wir verzichten hier auf eine allgemeine Diskussion des Periodenproblems der Ex- 
ponentialfunktion und des damit zusammenhängenden Mehrdeutigkeitsproblems für 
den Logarithmus. Im Hinblick auf das Folgende bemerken wir nur — mit den Bezeich- 
nungen von oben —: wenn p; ein Skalar von der Ordnung e; ist, dann wird 


Di +0 
(8, 2) I = 2 Asse 


eine Periode für e* sein, d.h. et = e#, Z und IT sind nämlich permutabel, und für 
zwei permutable Matrizen X und Y gilt nach $7, II die Beziehung 


(8, 3) eitY — eK .er —eF.eX, 
außerdem ist e!’!=1. Hieraus folgt, daß die Gleichung 
(8, 4) em A 
u.a. die Lösung hat 
(8, 5) X=ArllogP + IA. 


Wenn & eine komplexe Variable und A eine konstante Matrix ist, sieht man, daß 
die Funktionen e*4 und x = e4!%= die Differentialgleichungen 

(8,6) Y’(x)=AY(z) = Y(x2)A, bzw. zY’(z) = AY(z) = Y(x)A 
befriedigen ?). 

Da wir besonders x noch häufig verwenden werden, geben wir dessen kanonische 
Formen an. Ist die Normalform von A die folgende 


PR---0 
(8, 7) Para +0-( PRTRPENR ') 
ce 
wobei P; = o; + ©,, so hat man?) 
e—1 v 
u ee (log x) „» 
AA=ı' = P | 5 [@) 


Ey 
(8, 8) -( ) 
Erz 


Ist weiterhin 





| i 
= m ı A, 
(8, 9) ((— A)" = = 
2% (e — 0;) e 
dann wird 
wi (log x)” 
A — ee Ka A i 
(8, 10) x =2 4 we. 


Kapitel 11. 
Das Produktintegral. 
$ 9. Die Existenz des Produktintegrals. 


In der komplexen Ebene sei ein reguläres Kurvenstück £ ohne Doppelpunkte 
gegeben, und ?(x) sei eine Matrixfunktion, die längs dieser Kurve stetig ist. Wir suchen 





1) Peano [1]. 
2) Volterra [2], S. 30—31. 
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dann eine Matrixfunktion Y(x), die längs £ differentiierbar ist und die Differential- 
gleichung 


(9, 1) Y'(2) = Y(z) P(«) 
befriedigt. 

Bezeichnen a und x zwei beliebige Punkte auf £ und sind a, = a, @,,.. -, Au-ı, Am = X 
Punkte zwischen a und x auf der Kurve, so ist infolge (9, 1) 

(9, 2) Y (a„+1) — Y(a,.) = (Y(a,) P(a,) + H(a,, ö,)) 64; 
indem wir ö, = a„+1 — a, und 


Ya +6)— Y(&) _ 


H(2, 6) _ Ö 


Y'(x) 





setzen. 
Mit Hilfe der Gleichung (9, 2), die man auch 
Y(a.+1) = Y(a,) (1 + Pla.) 6.) + H(a,, 6.) 6, 
schreiben kann, läßt sich Y(x) durch Y(a) und die Größen H(a,, 6.) ausdrücken; eine 
leichte Rechnung zeigt, daß 


m—1l m—1 


m—1 
(9,3) Y(z) = Y(a) I (1 + P(a,) 6.) rs H(a,, 0,) 6, 11 (1 + P(a)ö,). 
Hieraus kann man schließen: falls eine längs & stetige und differentiierbare Lösung von 
(9, 1) existiert, wird sie dargestellt durch den Grenzwert 


9,4) Y(z) = Ya) -lim TZ (1 + P(a,)ö,), 


m>nr= 
wo das größte |ö„| > 0, wenn m «. 
Erstens ist nämlich Y’(x) nach (9, 1) auf £ stetig und somit gleichmäßig stetig, 
so daß 


ö 
Hz, 9)= 3 | (r@+9— Yaydı=o, wenn ö>0!). 
0 
Zu jedem hinreichend kleinen ö können wir demnach ein e bestimmen, so daß 
(9, 5) IH(z,6)| <eE 
für jedes x auf £ und lime =. 
dö0 
Zweitens können wir einen Skalar c so wählen, daß 
|P(a)| < eE 


für jedes x auf £, und hieraus folgt (vgl. (7, 15)) 
m—1 


m—1l 
(9,6) ET 1 + P(ia)Ö)|SE-M(i+cjö|E)seM.E<e*.E, 
u 


wobei / die Länge des Kurvenstückes von a bis x bezeichnet. 
Aus den zwei Ungleichungen (9, 5) und (9, 6) folgt 


m—1 


IY(x) — Y(a)II (1 + Pla,)d,)| Se E-0, wenn m> ©, 
und damit ist (9, 4) bewiesen. 


Daß andererseits der Grenzwert in (9, 4) existiert, wenn ?(x) stetig ist, kann man 
beweisen auf Grund der Identität 


m— 


1 m 
IT (1 + P(a,)d,)=ZTW”, 
I, = ZPia)... Pi) du... I =1, 


u? 


(9, 7) 





1) > bezeichnet: strebt gleichmäßig gegen. 
10* 
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wo die Summation sich über 


(9, 8) Ss, <<. <<, sm—i 
erstreckt. Es gilt nämlich für jedes w 
(9, 9) lim 77° = [ Pit)... Plt,) di,...di„= T,(a,2), 


wo der Integrationsbereich durch die Ungleichungen 

(9, 10) ash, SWS... - SUuST 
angegeben wird. Mit der Ungleichung 5 < c meinen wir dabei, daß beim Durchlaufen 
von £ in der Richtung von a nach x, b vor c passiert wird. Weiterhin ist 


ITW|S#-Z|6,...6,|-E, 


und hieraus folgt 


7} 
(9,14) im CD .E. 
u! 
Entwickelt man (|ö,| + + |öm-ı|)” nach der Polynomialformel, so werden die | | 
Koeffizienten von |ö,...ö,,] nämlich w!, so daß 
(9, 12) El... 8,18 8014 + ll), 
welches zusammen mit &|6,] <2 (9, 11) ergibt. Läßt man m — © streben, so bekommt i 
man schließlich ; Ä 
u 
(9, 13) Tas -E. 
Hiernach kann man leicht zeigen, daß 
m—1 ) | 
(9, 14) lim IT (1 + P(a,)6,) = T,(a,). 
Diesen Grenzwert, das Produktintegral von 1 -- P(t)dt von a bis x längs &, bezeichnen | 
wir mit 
N 
(9, 15) Y(a,2) = / (1 + Pit) dt). t 
Dieser wurde zuerst von Volterra [1, 2, 3] und später von Schlesinger [1, 2, 3] betrachtet, | l 


während die Reihenentwicklung (9, 14) für eine Lösung der Differentialgleichung (9, 1) 
von Peano [1] herrührt und später von Baker [1] wiedergefunden wurde. Eine allgemeinere 
Reihe wurde von Schlesinger [1] S. 41 angegeben, aber sie geht aus der von Peano durch 
eine Transformation hervor (vgl. $ 19). 

Der hier geführte Existenzbeweis für das Produktintegral gilt offenbar auch, u 
wenn P(x) nach Riemann integrierbar ist, und er kann leicht dahin ergänzt werden, daß 
man nur voraussetzt, daß P(z) beschränkt und nach Lebesgue integrierbar ist. Da wir 
aber im folgenden keinen Gebrauch davon machen werden, begnügen wir uns damit, auf 
die diesbezügliche Arbeit von Schlesinger [4] hinzuweisen. 


— 


$ 10. Die Rechenregeln des Produktintegrals. 


Aus der Definition des Produktintegrals folgt unmittelbar, daß, wenn a, b und c auf | d 
£ liegen, die folgende Relation besteht 


b e . 


(10, 1) ik + P(t)dt) fa + Pt)dt) - [u + Pit)dt). Ri 
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Setzen wir 


(10, 2) fa + Pli)di) =A, 


dann gilt diese Beziehung auch für c= a, und hieraus folgt 


E Mc 
(10, 3) (/ + Pros) = / (1 + Pt)dt). 
a b 
Liegen die Punkte a, zund x + hauf £, so ist infolge (10, 1) 
Y(a,c-+h) = Y(a,2) Y(z,x-+h), 
und infolge (9, 15) hat man 


z+h 
Mer) a [ Po— Pina IT T, 2 





Das erste Glied auf der rechten Seite geht offenbar gegen 0 mit A —0, und dasselbe gilt 
auch von dem zweiten Glied, da es nach (9, 13) absolut genommen 


ist, wobei / jetzt die Bogenlänge von x bis x -1- h bezeichnet, also kleiner ist als eine Kon- 
stante -h. Hieraus folgt 

. Y(a,2z+h)—!la,2) _ oY(a, z) 2) _ 
che u h a Yu 





= Y(a,x) P(x), 


und damit haben wir gezeigt, daß Y(a, x) als Funktion von x differentiierbar längs £ ist 
und die Gleichung (9, 1) befriedigt. 

In analoger Weise kann man leicht zeigen, daß Y(a, x) auch als Funktion von a 
längs £ differentiierbar ist, und die Gleichung 


(10, 5) oY(a, 2%) 


ca 
befriedigt. 
Fassen wir diese beiden Resultate zusammen, so können wir sagen: Das Produkt- 
integral als Funktion der zwei Veränderlichen a und x besitzt das totale Differential 


(10, 6) dY(a, x) = — P(a) Y(a, x) da + Y(a, x) P(x) dx. 


= — P(a) Y(a, x) 


$ 11. Die Peanosche Reihe. 


Wir werden jetzt die Peanosche Entwicklung des Produktintegrals etwas näher 
untersuchen. Man sieht leicht, daß das u-fache Integral in (9, 10) als zwei iterierte 
Integrale 


(14,1) Ta, 2) =[ SS Po): Pl)du---di, 


at ee 


z iu 


u ca f Piu) ++ Pl)dty +: di, 
dargestellt werden kann, und hieraus liest man die zwei Rekursionsformeln !) ab: 


(11, 2) T +1la, 2) = f P(t)T ‚(t, z)dt = f T ‚(a, t) Plt)dt. 





1) Peano [1]. 
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Wird das erste Integral nach a und das zweite nach x differentiiert, so findet man 
die partiellen Ableitungen von 7,+1(a, x) und hieraus wieder das totale Differential 


(11, 3) AT „+ıla, x) = — P(a)T,(a, x)da + T,(a, z)P(x)dx. 


Wird die erste Formel in (11, 2) mehrere Male nacheinander angewandt, dann be- 
kommt man 


T u+,+1(a, x) = f vr [P) UM Pitu+ı) Toltu+ı, 2) dty+ı PO dt, 
a t, 


und durch Vertauschen der Integrationsordnungen 


7 Eu+1 t: 


T u+,+1(@, x) = [ S Pr -[ Pit) TOR Pltu+ı) Toltur, x) dt, di: diu+1; 
woraus folgt 
(11, 4) T u+»+1(@, x) = [T,(a, t) PT Ct, x) dt. 


Mit Hilfe dieser Formel kann man einen einfachen Ausdruck für das Restglied der 
Peanoschen Reihe finden. Setzt man 


(11,5) Y(a, 2)=&T,(a, 2) + Rula, 2), 
so hat man nämlich 
Rula, 2) = z ER Vo f B> T (a,t) Pit) T„(t, 2)dt, 
und hieraus folgt 
(11, 6) Rn(a, x) (Ye, t)P(t)T(t, z)dt. 
Analog hierzu ist !) 


(11,7) Rn(a, x) - (T,(a, t) P(t)Y(t, z)dt. 


$ 12. Die inhomogene Gleichung. 
Die inhomogene Gleichung 
(12,1) Y'(z) = Y(a)P(z) + R(e), 
wo P(x) und A(z) auf £ stetig sind, ist leicht zu lösen, wenn man ein Fundamental- 
integral zu (9, 1) kennt. Man hat nämlich in Analogie zu (9, 3) 


m—1 m—1 m—l 


Y(a) = Ya) IL (1 + Pla) + 2,(Rla,) + H(an 64))84 IL (1 + Pla), 


so daß die Lösung, die O ist, für x = a, dargestellt wird durch 


m—1 m—l 
(12, 2) Yı(a, x) = lim 3, R(a,)ö, IL (1 + P(a)ö.). 


Da nun das Produkt (9, 4) gleichmäßig gegen das entsprechende Produktintegral kon- 
vergiert, so wird 


1 


h7 (1 + P(a,)ö,) —/ (1 + Plt)d)=0, wenn m>o», 
u Fü 


!) Vgl. Schlesinger [1], S. 41, Formel IV. 








al 


di 
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und daraus folgt 
(12, 3) Yı(a, x) = [ R(t)Y(t, a)dt. 


Analog hierzu wird die Lösung von 
(12, 4) Y'(x) = — P(x)Y(x) + Ki), 
die O ist, für x = a, dargestellt durch !) 


(12, 5) Y,(a, x) (Ya, Rlı)dt. 


Setzt man z. B. in (12, 3) die Peanosche Reihe mit dem Restglied (11, 6) für Y(t, x) 
ein, so findet man 


m 


(12, 6) Y,(a, x) = S,(a, x) + Ru(a, x) ’ 
wobei 


(12,7) S,(a, x) = ROT,dı, x) dt 
und 
Rn(a, x) nd Ru) Y(t, r)P(r)T „(t,x)drdt 
(12, 8) | = 
—= [ Yı(a, z)P(r)T (rt, z)dr, 





\ 


oder bei Anwendung der Formel (11,7) für R,„ 


(12, 9) Rla, x) = (Smta, t) P(r)Y(tT,x)dr. 
Schließlich leitet man leicht aus 11,4) die Formel 

(12, 10) Su+,+1(a, X) = [5,(a, t) P(t)T,(t, x) dt 
ab, die auch für «= —-1 gilt, wenn wi S_ıla, 2) = R(x) setzt. 


Für »=0 ergibt sich die Rekursionsformel 
(12, 11) Su+1(a, x) = [S,(a, t) Pit)dt, 


d. h. die Folge (S,(a, x)) ist durch genau dieselbe Rekursionsformel bestimmt, wie (T,.(a, x)), 


nur lautet das nullte Glied der zwei Folgen S,(a, x) = f R(t) dt bzw. T,(a, x) =1. 


Zur Anwendung der Formel (12,3) bemerken wir: Sind Y(a, x) und Z(a, x) die 
Produktintegrale von 1 + P(t)dt bzw. 1 + Q(t)dt, dann ist Y(a, x) — Z(a, x) die Lösung 
der Gleichung 


(12,12) V’(2) = V(2)Q(2) + Y(a, x) R(x), P(x) = Q(x2) + R(x), 
die O ist, für c=a. Folglich ist 


| Y(a, 2) — Z(a, x) = f Y(a, t) R(t) Zit, x) dt 
(12, 13) | " 





n f Z(a, t) R(t) Y(t, x) dt. 


a) Peano [1]. 
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Ist P(x) weiterhin von einem Parameter A abhängig und hat es eine in x und } 
stetige partielle Ableitung nach t, so kann man hieraus leicht einen Ausdruck für den nach 
A partiellen Differentialquotienten des Produktintegrals 





N 
(12, 14) Y(a,x; A) = | (1 + Pit: A) dt) 
ableiten. Man findet 
OYla, 2; 2) _ [ ndPl;A) ,,, 
(12, 15) Zu Y(a, 2; 3). Yu, a; A)dt. 


In dem speziellen Fall 
(12, 16) P(xz; 4) = AP(x) 
wird die Formel besonders einfach: 
(12, 17) al E -[ Y(a,t; }) P(t)Y(t, x; A) dt. 


4 


$ 13. Fundamental- und Partikularlösungen. 


Eine Lösung von (9,1) nennt man eine Fundamentallösung, wenn jede andere 
Lösung ein Produkt dieser Lösung und einer konstanten Matrix ist. 
Da man nach (9,4) schreiben kann 


(13, 1) Y(z) = Y(a)Y(a, 2) 
und dies für jede Lösung gilt, ist Y(a, x) eine Fundamentallösung. 
Die Bedingung dafür, daß Y,(x) eine Fundamentallösung ist, kann auf verschiedene 
Weise ausgedrückt werden. Einmal kann man eine Konstante A finden, für die 
Y(a,z) = AYo(2); 
aber diese Konstante muß infolge (13,1) gleich Yo (a) sein, das also existieren muß. 
Existiert Y5 (a), dann ist infolge (13, 1) 


(13, 2) Y(a, x) = Yo'(a)Yo(2) 
und deswegen für eine beliebige Lösung 
(13, 3) Y(z) = Y(a)Yo'(a)Yo(2). 


Es ist also die Existenz von Yo ‘(a) notwendig und hinreichend dafür, daß Y,(x) eine Fun- 
damentallösung ist. Da indessen a eine beliebiger Punkt auf £ sein kann, so läßt sich 
die Bedingung auch so ausdrücken: es muß Y, (x) existieren, d. h. Det. Y,(x) # 0 sein 
für jedes x auf £, oder auch: es darf keine Konstante A # 0 geben, für die AY,(x) = 0 
ist, in einem Punkte auf £. 

Die Existenz von Y”’(a, x) geht aus (10, 3) hervor, und daß Det. Y(a, x) + 0 ist, 
folgt selbstverständlich hieraus, aber kann auch durch direkte Ausrechnung der Deter- 
minante bewiesen werden. Es ist nämlich 


Det. Y(a, 2) = lim IT Det, (1 + P(a,) 8.) 
und 
Det. (1 + Pt)d)=1 +2 pP; 6 + + Det. Pl) =", 
wobei n(t,6)>0, wenn ö —0; folglich ist 


(13, 4) Det. Y(a, x) = exp. j£ p,(t) dt) , 


also überall dort #0, wo P(x) stetig ist. 
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Diese, von Jacobi [1] angegebene Formel, wurde beträchtlich verallgemeinert von 
Lappo Danilewski [1]. 

Ist Y(x) nicht invertierbar und ev. rechteckig, dann hat sie infolge (13, 1) den- 
selben Rang r für alle x, für die & p,,(x) stetig ist. Wir nennen sie dann eine Partikular- 
lösung vom Range r. Durch Multiplikation von links mit einer passenden Vertauschungs- 


matrix wird sie auf die Form ei gebracht, wo Y, nach rechts invertierbar ist und die 
2 


Zeilen in Y, aus den Zeilen von Y, abgeleitet werden können. Y kann also durch eine 
Lösung von der Ordnung (r,n) ersetzt werden. Diese kann weiterhin in der Form 
(1,, On-r,r) : Yo(x) geschrieben werden, wobei 1, die Einheitsmatrix von der Ordnung r 
und 0,-,,r die Nullmatrix von der Ordnung (n — r,r), während Y,(x) eine Fundamental- 
lösung ist. Durch Multiplikation von rechts mit einer passenden Vertauschungsmatrix 
läßt sich nämlich Y,(a) auf die Form (C,,, Cj,) bringen, wo C,, invertierbar ist, und 
demnach zerlegt werden kann in 


Ca “n) 
TE ya 


wobei der zweite Faktor invertierbar ist. Das Produkt dieser, der Vertauschungsmatrix 
und Y(a, x), ist somit eine Fundamentallösung. 
Für eine Lösung von der Ordnung (1, n) gebrauchen wir zuweilen die Bezeichnung 


Vektorlösung. 


$ 14. Die Konvergenz des Produktintegrals. 
Ist P(x) auf £, ausgenommen im Punkte a, stetig, so nennt man das Produktintegral 


von a bis x konvergent, wenn 
z 


(14,1) Ju + Paya)=tim [1 + Poyaı) 
u nz 


existiert und invertierbar ist. Die Bedingung dafür, daß dies der Fall ist, lautet: 
4 
(14, 2) [a + Pway—1 <eE 
8 | 
| $ı I 


für alle &, und £, zwischen a und x, die hinreichend nahe bei a liegen. 
Daß diese Bedingung notwendig ist, folgt aus (5,5) und aus 


Y (En &) —1 = (VE 2) — Fl 2)) 1, $2)- 
Daß sie auch hinreicht, erkennt man aus den Relationen 
Y(&, ) un Y(£,, z) u (FE &) Er 1) Y(&, z) 
und 
Y(z, &) — Ya, 5) = Ya, E29) Ye) —2). 
Diese zeigen einmal, daß sowohl Y(£, x) als auch Y(x,£&) für £—>a beschränkt sind, 
dann, daß Y(£,, 2) — Y(&,, x) und Y(x, &) — Y(z, &,) beide — 0, wenn £, und £, unab- 
hängig voneinander gegen a konvergieren. 
In der Praxis reicht oft das folgende Kriterium aus: 


I. Ist [ P(t)dt absolut konvergent, so ist auch das Produktintegral konvergent und wird 


a 
durch die Peanosche Entwicklung dargestellt !). 
1) Volterra [1], S. 54—55 beweist die Existenz des Produktintegrals unter der Voraussetzung 


cE 
EOS Tan ö>0. 
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Die Peanoschen Entwicklungen für Y(x,&) und Y(£, x) sind nämlich gleichmäßig 
konvergent, da sie offenbar absolut konvergente von & unabhängige Majorantenreihen 
besitzen. 


$ 15. Das unbestimmte Produktintegral. 


Eine gegebene Lösung von (9,1) kann natürlich jederzeit dargestellt werden als 
Y(a) - Y(a, x), mit passend gewähltem a, aber in gewissen Fällen ist eine andere Form 
der Darstellung wünschenswert. 

Ist die Folge (T,(x)) eine beliebige Lösung des Rekursionssystems 


(15,1) T,(z) = T,„-ı(z)P(x), T,(z) konstant, 
so wird die Reihe 
(15, 2) Y(z) =2 T,(«) 


augenscheinlich eine Lösung von (9, 1) darstellen, wenn sie gleichmäßig konvergiert. Das 
ist aber der Fall, wenn Konvergenz in einem einzigen Punkt stattfindet. Durch suk- 
zessive Integrationen von (15,1) für « = 0,1,2,... findet man nämlich 


(15, 3) T,(2) =2 T,(a)T,-,(a, 2), 


so daß (15,2) als das Cauchysche Produkt von Y(a)=&T,(a) und Y(a, 2) = £T,(a, 2) 


aufgefaßt werden kann. — Als Spezialfall von (15, 3) ergibt sich offenbar eine neue Re- 
kursionsformel zwischen den Funktionen T,(a, x): 


(15, 4) T,.(a, 2) =& (a, b)T,-,(b, 2). 


Inwieweit das so gebildete „unbestimmte Produktintegral‘‘ eine Fundamental- 
lösung ist, muß in jedem einzelnen Fall für sich untersucht werden. Doch für eine gewisse 
Klasse von solchen Bildungen läßt sich diese Frage allgemein entscheiden. Bildet man 
die 7,(x) = (&£(x)) durch die sukzessiven Integrationen 


(15, 5) (x) = J Zur) p„(t) dr 
mit denselben endlichen unteren aa für jedes «, so hat man 
t, 
( 
(15, 6) ik (2) = 2 RA en Jan - [ p,(t) P;,,;,(t,) “wr P;,_..(,) dt, . 


Durch eine Schlußfolgerung analog zum Existenzbeweis in $ 9, erkennt man nun, daß das 
hierdurch definierte Produktintegral als der Grenzwert 


m—1 
(15, 7) lim II A+ (p„(ar)ö@)), oe — ai alt _ ee 
dargestellt werden kann. Dabei sind a, u =1,...m—.1 Punkte auf einem Kurven- 
stück, das a) = a, mit a) = x verbindet. Indem wir A = (a,,) setzen, schreiben wir 
für das besagte Produktintegral 


(15, 8) Yaz) = | (1 + Pü)dı) = | (1 + (put) 
(4) (4) 


und weiterhin 


(15, 9) 7,02) = [ Tui) Pi) de. 


(4A) 








le 
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Analog dem Beweis für (13, 4) beweist man leicht die entsprechende Formel 


(15, 10) Det. Yı(z) = exp. (ES put), 
woraus die Invertierbarkeit von Yy,(x) ersichtlich ist. — Übrigens kann man auch 
leicht Yz7'(z) bilden: 
0 
(15, 11) Ya (2) = lim II (1 — (pu(aiw) 8%), 
und wir schreiben hierfür 
(4) 
mM 
(15, 12) Yz!(z) = (1 — Pit) dt). 


Die Glieder seiner Peanoschen Entwicklung sind bestimmt durch die Rekursions- 


formel 
(4) 


(15, 13) T,(z) = Pt) T,-ı(t) dt. 


Für die Diskussion eines asymptotischen Problems der skalaren Gleichungen n-ter 
Ordnung verwendet Horn [3] S. 175 eine sukzessive Approximation, die mit der Bildung 
az0, ısk 


o,i>k gewählt wird. Man vergleiche 


von Yı(x) äquivalent ist, wenn ai: -| 
übrigens auch $ 27. 


$ 16. Die adjungierte Gleichung. 


Im folgenden soll die Variable auf dem oben besprochenen Kurvenstück & liegen, 
und alle vorkommenden Funktionen sollen der Einfachheit halber, als überall differentiier- 
bar angenommen werden. Weiterhin sollen, um die Formeln so übersichtlich wie möglich 
zu machen, dort, wo Mißverständnisse ausgeschlossen sind, die Variablen weggelassen 
werden, so daß im allgemeinen z.B. P und Y statt P(x) bzw. Y(x) oder P(t) bzw. Y(t) 
geschrieben wird. 

Aus der Identität 


(16, 1) (ST' =S’T+ST 

folgt für S= T", 
(16, 2) 1) r’r m”, 

und hieraus erkennt man: wenn Y ein Fundamentallösung von 
(16, 3) !’=rYP 

ist, dann wird 
(16, 4) Tamm“, 

Die Transposition dieser Gleichung zeigt weiter, daß 
(16, 5) Z=—ZP=—zZxPı), 

die die zu (16, 3) adjungierte Gleichung genannt wird, die allgemeine Lösung 
(16, 6) ER 


hat, wobei C konstant ist. Diese Beziehung läßt sich auch folgendermaßen ausdrücken: 
YZ=-YxZ=(C, 


1) In diesem Paragraphen soll allgemein ein Stern bei einer Matrix deren transponierte Matrix bedeuten, und 
X die zeilenweise Multiplikation bezeichnen, 





11* 
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d. h. das zeilenweise gebildete Produkt von Y und Z ist eine Konstante. Speziell gilt für 
die Produktintegrale 
Y(a, 2) x Z(a,2) =1 


oder auch 
z 


(16,7) Y*(z,a) = Z(a, x) = Ta- P* dt). 


Wenn Y,Z und P beliebige Matrizen und W = YZ ist, kann man (16,1) auch 
folgendermaßen schreiben: 


(16, 8) W’ = (Y'’— YP)Z -- Y(Z’ + PZ). 
Ist speziell Z ein Integral der Gleichung 
Z=—PZ, 


dann bekommt man 
W' = (Y’'’— YP)Z. 
Bezeichnet außerdem Y ein Integral der inhomogenen Gleichung 
Y’'’=YP+R, 

so findet man wieder Formel (12, 3). 

Die Relation (16, 8) schreibt man gewöhnlich in der Form!) 

(16, 9) (Yx M)’ = (Y’— YP)xM-+TY’x (M-—- MxP) 
und stellt sie der Lagrangeschen Identität für eine skalare Gleichung r-ter Ordnung zur 
Seite. Befriedigt M die zu (16,3) adjungierte Gleichung (16,5), so bekommt man 

(Yx M)’=(Y'— YP)xM 

und nennt M einen Multiplikator zur Differentialform Y’— YP. 


$ 17. Derivation. 
Für die Derivation ?) 


(17,4) D,T= T"T 
gelten infolge (16,1) und (16,2) die Rechenregeln 

(17, 2) D«ST)= T"'-D,$S:-T+D,T 
und 

(17, 3) DT'"=—- TT'"=-—-T.DT.T". 
Speziell wird 

(17, 4) D,(AT)=D,;,T, 


wenn A eine Konstante ist. 
Die Bedeutung dieser Operation beruht darauf, daß sie die Umkehrung der Produkt- 
integration ist; ist Y ein Fundamentalintegral zu (16, 3), dann wird nämlich 


(17,5) D,Y=P., 
so daß 
(17,6) Y'(a) Y(z) = Ya 2) = fa -—- D, Ylt) dt). 


$ 18. Die Q-Transformation. 
Durch die Transformation 


(18,1) Y=ZT, Z=YT" 


E Vgl. Schlesinger [1] S. 33. 
*) Volterra [1], 5.1734 und Schlesinger [1], S. 26—28. 
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geht die Gleichung (17, 5) infolge (17,2 Jund (17, 3) über in 


(18, 2) Du P*: 
für den Ausdruck auf der rechten Seite führen wir das Symbol 
(18, 3) APIT)=TPT"— TT"= T(P—D,T)T" 


ein und nennen ihn die Q-transformierte von P durch T. Demzufolge kann man (18, 2) so 
schreiben: 


(18, 4) Z’=ZaXP|T). 
Für diese Transformationen gelten einige einfache Rechenregeln. 
So ist 
(18, 5) Q(0|7)= D,T" 
und 
(18, 6) AP|A)=APAT, 


wenn A konstant ist. Die Q-Transformation verallgemeinert also die Derivation und ın 
gewissem Sinne auch die Transformation (2,1), die dem Äquivalenzbegriff!) zugrunde liegt. 
Wird in (18, 4) 


Z = US 
gesetzt, so hat man im ganzen 
Y=U-.ST 

gesetzt und findet die Kompositionsregel 

(18,7) aaXP|T)IS)=RWP|ST). 

Besonders erhält man für $ = 7” ein gewisses Reziprozitätstheorem: Aus 

(18, 8) APIT)=0 folgt XE|T")=P. 
Es geht dies im übrigen auch aus der Beziehung 

(18, 9) aPpıT)=T"PT+-T"T 
hervor. 

Schließlich sei bemerkt, daß 

(18, 10) AXP-+Q|T)=XP|T) + TOT", 
und 

(18, 11) D.(ST) = XD,S | T) 
ist, sowie, daß (16, 3) auch 

(18, 12) QaP|Y)=0 


geschrieben werden kann. 


$ 19. Transformation des Produktintegrals. 
Ist Y eine Fundamentallösung von (16, 3), so hat man infolge (18, 1) 
Y*'(a) Y(z) = T'(a)Z”'(a) Z(x) T(x) 
oder, bei Einführung des Produktintegrals ?) 


T zT 


(19, 4) fa + Pdt) = 7a) [di + 2P| T)dt) T(e). 


a 





1) A. Loewy [1] hat eine tiefgehende Untersuchung der Äquivalenz vorgenommen, nach der zwei Matrizen ? 
und Q äquivalent sind, wenn die eine aus der anderen durch eine Q-Transformation hervorgeht, wobei die Elemente 
von T in demselben Rationalitätsbere'ch liegen wie die von P. 

2) Vgl. Volterra [1], S. 50. 
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Eine andere Form bekommt man dadurch, daß man P durch P +0 ersetzt und 
für T eine Fundamentallösung von (16, 3) wählt. Durch Anwendung der Additions- 
formel (18, 10) und von (18, 12) erhält man 


T E7 


19,2) [a +1(P+Qa) = T*%a)| (+ TOT'a) Te). 


Setzt man in die rechte Seite die Peanosche Reihe für das Produktintegral ein, so 
ergibt sich die Reihenentwicklung für eine Lösung der Gleichung 


(19, 3) !"=Y(P+Q), 
die Schlesinger [1], S. 36—45, durch die Methode der sukzessiven Approximationen ab- 
geleitet hat. Dies geschieht so, daß man 


(19, 4) Yo z Y, 
setzt und die Glieder der Reihe durch die Rekursionsformel 
(19, 5) Hn=HP, = T,P+T7Y.ı0 


bestimmt, indem die inhomogenen Gleichungen durch die Formel (12, 3) mit gleichem a 
für alle w gelöst werden. 


Wir werden gelegentlich die Reihe in einer anderen Form benutzen, indem wir die. 


Funktionen 
(19,6) Vi) = | + 7, QM) Ta, nd), Ta, = T"TW, 


einführen, so daß 

(19,7) Y(a,x) = T(a,x)V(x). 
Durch Anwendung der zweiten der Rekursionsformeln (11, 2) für die Glieder in der be- 
sagten Peanoschen Reihe findet man 


(19, 8) Vea)=£ Tue), 


wobei 7,(2) = 1 und 


(19, 9) T „+1(&) =/[ T(x, t) T (t) Q(t) Tt, x) dt ’ u > 0. 
Geht man von der differentiierten Rekursionsformel aus, so bekommt man 
(19, 10) Tı+ı — Tarı P + PTırı = T,Q. 
Dem entspricht, daß V die Gleichung 
(19, 11) V’’=V(P+Q)— PV 
befriedigt. 


8 20. Die Transformationsgleichung. 
Will man eine gegebene Gleichung (16, 3) in eine andere von derselben Form 


(20, 1) Z=ZQ 

transformieren, so handelt es sich offenbar darum, eine Lösung von 
(20, 2) AıPpıV)=Q 

oder 
(20, 3) V’’= VP—QV 


zu bestimmen. Gleichungen dieser Art spielen deshalb bei vielen Problemen eine be- 
deutende Rolle. 
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Durch die Transformation 


(20, 4) V=S'WT 
geht (20,3) über in 
(20, 5) Ww' = WAP | T)— XQ | S)W. 


Es ist nämlich 
QAXPIWT)=NMP|SV), 
und nach der Kompositionsregel (18, 7) ist die linke Seite 
= AXP|T)W), 

während die rechte Seite 

= AP V)IS)= UQI|S) 
ist; aus der Identität beider Ausdrücke folgt (20,5) !). 

Als Spezialfall dieser Transformation sei hervorgehoben, daß 


(20, 6) Y = 72 
die Gleichung (16, 3) überführt in 
(20, 7) Z=ZP-—D.T-Z. 


Mit Hilfe von (20,5) kann man die Form der allgemeinen Lösung von (20, 3) be- 
stimmen. Bedeuten nämlich 7 = Y und 5 = Z Fundamentallösungen der Gleichungen 
(16,3) und (20,1), dann nimmt (20,5) ie Form 

W=0, d.h.W=tC 
an, so daß die allgemeine Lösung 

(20, 8) v=ZT’cy 
ist, wobei C eine willkürliche Konstante bedeutet ?). 

Für eine gegebene Lösung kann man die Konstante dadurch bestimmen, daß man 
x = a setzt; es ergibt sich dann 

(20, 9) V(x) = Z(z, a)V(a)Y(a, x). 

Besonders sei bemerkt, daß die allgemeine Lösung der Gleichung 

(20, 10) U'=UP—PU, 
mit der wir uns weiter unten noch beschäftigen werden, 

(20, 41) U= Y’CY= Y(z, a)U(a)Y(a, x) 
ist; jede Lösung von (20, 10) ist also einer Konstanten äquivalent. 


Ist V eine gegebene Lösung von (20, 3), dann kann man nach dem vorhergehenden 
5, T und C so bestimmen, daß 


(20, 12) V=SCT. 
Infolge (20, 5) wird 
(20, 13) .  CAP|IT)= RE|S)C, 


und weiterhin ist leicht zu erkennen, daß, wenn S, T und C diese Relation befriedigen, 
V eine Lösung von (20, 3) ist. Man braucht nur die Gleichung 


0=(C’"=CRP|T)— RXQ|S)C 
dadurch zu transformieren, daß man 
C=sSVT" 
1) Es sei angemerkt, daß J,agranges Identität (16, 8) auch so geschrieben werden kann: 
Ww=waPpız')—aPpıyWw, W=YZ 


und somit als ein Spezialfall von (20, 5) erscheint, nämlich für@= Pund V=1. 
2) Schlesinger [2], S. 291. 
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setzt; man kommt dann infolge (20,5) zu (20,3). (20,12), ergänzt durch die Bedingung 
(20, 13), die sich für $ =Z, T = Y auf eine Identität reduziert, ist also eine vollständige 
Lösung der betrachteten Gleichung in allgemeinerer Form als (20, 8). 

Es mag noch bemerkt werden, daß, wenn 


(20, 14) rT=07, ZW 
gesetzt wird, ® und Y durch die Relation 
(20, 15) cd = YC 


verbunden sind. 


Ersetzt man in Gleichung (20,3) Q durch @— RV”', dann geht sie über in die 
inhomogene Gleichung 


(20, 16) V’’=VP—QV-R. 
Diese läßt sich durch die Transformation (20, 4) überführen in 
W' = Wa(P| T)— 2(@ — RV|S)W, 
die aber unter Anwendung der Additionsformel (18, 10) zu 
(20, 17) W' = Wa(P| T)— 20 | S)W -- SRT | 
umgeschrieben werden kann. Wählt man ebenso wie in (20,8) 7=Y,5=Z, so be- \ 
kommt man 


W'=ZRY", 
so daß die inhomogene Gleichung (20, 16) die Partikularlösung 


ET RE da 





(20, 18) | V.{a,2) = [Z, t) Rit) Y(t, z)dt | 


hat. Eine beliebige Lösung kann infolge (20, 9) 
(20, 19) V(z) = V,(a, x) + Z”'(a, x)V(a) Y(a, 2) 
geschrieben werden. (20, 18) enthält die beiden Formeln (12, 3) und (12,5). 


$ 21. Das Reduktionsproblem. | 
Die Integration einer Differentialgleichung kann man als ein doppeltes Problem & 
auffassen: erstens die Transformation der Gleichung in die möglichst einfache Normal- # 


form, wobei die Transformatoren einer gegebenen Funktionsklasse angehören, und 
zweitens in eine darauf folgende Diskussion der Integration der Normalform. | 
Im allgemeinen geht das Reduktionsproblem darauf hinaus, eine Normalform 


P,(x2) 0 ---0 
(21,1) P(x) = | EEE ) 
0 0... P„(x) 
hervorzubringen, wobei die P;(x) nicht weiter reduzierbare quadratische Matrizen sind. 2 
Im günstigsten Fall sind sie alle von erster Ordnung, und man spricht dann von voll- 
ständiger Reduktion. 
Es wird in theoretischen Untersuchungen oft übersichtlicher sein, die Reduktion 
schrittweise auszuführen, indem man zuerst P auf die Form 





EEE EEE ETF ITETUTTTNG ne 


0 
24,2 3" Pr ) 
(21, 2) (6.0 
bringt und danach 0, und Q, auf dieselbe Art so weit wie möglich zu reduzieren versucht. & 


Qı 
0.Q 


i reduziert und multipliziert man diese mit dem Transformator, der P 


Reduziert man z.B.0Q, durch den Transformator T, zu | x ‚sowirdQ@ durch den Trans- 
1 


T, 


0 


formator ( 








a 
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in Q überführte, so haben wir eine Transformation, die P auf die Form 


0 00 
(21, 3) (0 Q, 0) 
00 
bringt. Und so kann man offenbar weitergehen, bis die Form P(x) erreicht ist. Das 


fundamentale Problem ist deshalb die Reduktion von P auf die Form 0, die relative Re- 
duktion. 


$ 22. Transformation mit Relativmatrizen. 


Im $1 wurde gezeigt, daß ein Transformator, der ja invertierbar ist, in das Produkt von 
einer Vertauschungsmatrix V und dreirelativen Elementarmatrizen zerlegt werden kann: 


A )(t 0)(f T,) 
.. r=|(, To ı )F 


wobei 7,, quadratisch von der beliebigen Ordnung r < rn ist. Im Hinblick auf die Kom- 
positionsregel müssen wir uns deshalb über die Wirkung einer Transformation mit jeder 
von diesen klar werden. 


Die Transformation mit V bewirkt nur die Vertauschung einiger Zeilen und 
Spalten in ? untereinander, aber so, daß die Diagonalelemente in der Diagonale bleiben. 


Indem man P in (5 5) zerlegt und beachtet, daß 
Pa Pa 


a (1 


bekommt man die wichtigen Transformationsformeln 
(22, 3) al? | ( =) . (* + Ta Pa Pa—PuTet Taf — TıaPaı Ta — a) 


'\04 er Pan Pa — Pa Tıa 

1 0 Pa—PaT P | 
22,4 ap ))= ( 11 12 m a - ; 
| Tzı 1 Pa—PaTa+ Ty Pı— Ta Pı2 Taı — Tai, Pa+ Pu) 


(22, 5) a(pi(t" 0 ))= bb In), Tıı Pam, 
'\0 T a2 T a9 Paı Tıı , QP| Te) 

Will man nun durch drei sukzessive Transformationen dieser Art eine relativ re- 
duzierte Matrix hervorbringen, so muß offenbar die rechte Seite in (22,3) relativ halb- 
reduziert sein, d.h. das Element in der ersten Zeile rechts muß Null sein. 7,, muß also 
die Gleichung 

(22, 6) Ta= Pat PuTıe — Tıa Pa — TıaPaı Tıa 
befriedigen !). Da die letzte Transformation die Reduktion nicht weiter befördert, muß, 
um die Reduktion zu vollenden, 7,, eine Gleichung von der Form 


(22, 7) T, = Ra —0Q2Ta + Tapıı 
befriedigen. 

Bei Anwendung dieses Prinzips muß man beachten, daß T,,, Tıa, Taı, Tag nicht 
immer zu derselben Funktionenklasse wie T zu gehören brauchen, aber wie in einer 
weiteren Arbeit gezeigt werden soll, wird man bei wichtigen asymptotischen Integrations- 
problemen 7, und 7,, in dieser Klasse bestimmen können, während 7,, und 7,, gleich 
1 gesetzt werden können. 





1) Vgl. Schlesinger [1], S.113, Gleichung (D). 
Journal für Mathematik. Bd. 171. Heft 2. 12 
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$ 23. Relative Reduktion vermittels Partikularlösungen. 


Kennt man zu Y’= YP eine Partikularlösung vom Range r, so kann man die 
Integration durch Lösung einer Gleichung von der Ordnung n —r und eine Quadratur 
ausführen. 

Ist Y,, eine invertierbare Teilmatrix von r-ter Ordnung, so transformieren wir 
die Gleichung mit 


_ (Ya, Yız) _ (Yu, O\f1, —. 
BE r-( 0,1 )-( 0, 6 1 
und finden nach (22,3) und (22,5) 
0,0 

23,2 AP|T)= -( E 

u a 
wobei 

(23, 3) Qa=PaYı, Q&= Pa— PaYı Yı. 
Setzt man 

A ; 

23, 4 Y=[( u u)= 27, 

Be u 
dann wird 

10 

23,5 Ei | 

aaa Zu Zu) 
also durch 

(23, 6) Zaı u Zee Op Zya Bas Zar O2 
bestimmt, und weiterhin ist 

(23, 7) Yaı en. Zaı Yın Y. Sn Ze 7 Zaı Y- 


$ 24. Produktintegration durch Quadratur. 
Wenn P/(z) und T,(a, x) permutabel sind, kann die Produktintegration auf die 
Quadratur zurückgeführt werden. Ist nämlich eine Funktion F(x) mit ihrem Differential- 
quotienten permutabel, dann wird 


ze ei — u(F(e))""'F'(x), 


und hieraus folgt durch Induktion 





1 oOT,(a,t 1 
241) Tl) = [Tail 5 dt = lTıla 2), 


so daß!) 


be) 4 r T [ruya 
24,2 Y(a,z)= 3 —(T,la, 2) Y=e"" = e 
(24, 2) (a, 2) = 3, ;Tıla, 2) 
Ein besonders wichtiger Spezialfall ist 
wobei A eine Konstante und p(x) ein Skalar ist. Dann ist 


z 


4 | p)d 
(24, 4) Y(a,z)=e* 


Für p(x) = 2 findet man 





1) Peano [1] hat die Formel (24,2) in dem wichtigen Spezialfall angegeben, wo P(t,) und P(t,) für be- 
liebige i, und t, auf dem Integrationsweg permutabel sind. 





ern ni 
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2 fü+.2)-(2)' 


Wird in (24,4) A = E gesetzt, so bekommt man unter Anwendung von (7, 15) eine für 
Majorantenbetrachtungen recht nützliche Formel 


N E fi pie) dt mo ar 
(24, 6) [u + pr Ed) = e® -1+(@ —A)E. 


Eine andere Möglichkeit für die Quadratur liegt vor, wenn sich ? als eine Relativ- 
matrix zweiter Ordnung von der Form 


pP (0 n 


schreiben läßt. Dann ist 
Pit) Pi) =, 


so daß die Peanosche Reihe bis auf die ersten zwei Glieder verschwindet. Man findet also 


£ A 2, 
(24, 7) J (i + A 0) d) -(; Pd, ) 


Ist ? relativ halbreduziert: 
(24, 8) Bi | 


und kennt man die Produktintegrale 


x 


24,9) Tu) = [+ Pudı), Tata, z)= [+ Padt), 


so kann man Ya, x) transformieren mit, 


und danach durch Quadratur bestimmen !): 


z 


Yla,2) = J (1 + kan 17 0) dı) Er 2) 


(24,10) | T ı(a, %) . | 
A [i Pan, d Palı) Tut, D)dt, Tala ) 





| 


Allgemeiner kann die Produktintegration durch Quadratur ausgeführt werden, 
wenn P eine Relativmatrix von der Ordnung m und der Form 


0.0 +-++-0 0 


ist. Dann verschwinden in der Peanoschen Reihe alle Glieder von dem m-ten an. Wenn 
P also relativ halbreduziert ist, d.h. von der Form 





1) Vgl. Schlesinger [1], S. 106—107. 
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Pı 0 ---0 
(24, 12) »-[P Pa : : +0 ) 


Pmı Pm2 ur Pmm 
so kann Y durch Quadratur bestimmt werden, wenn außerdem die Produktintegrale 


- 
T;(a, x) = + P;dt), a 


bekannt sind. Speziell kann die Produktintegration durch Quadratur ausgeführt werden, 
wenn P vollständig halbreduziert ist: 





Pıı 0 ... 0 
(24, 13) r(® Pr 0 ) 
Prı Pn2 Prn 
da ja 
P- fpict 
Ju+r,d)=e 
ist. 


S 25. Ungleichungen für das Produktintegral. 
Ist Q(t) eine für alle 2 auf £ reelle positive Matrix, für die 
(25, 1) ı P)]| SQ() 
ist, und setzt man 


Fa, 2) - Ta +owıaıı 





25,2 > 
= ZT,(a, x) + Rnla, 2), 
so hat man offenbar 
(25, 3) IT,(a,2)| < T,(a, x), u=0,1,2,... 
und damit 
(25, 4) IY(a,x)| < Y(a, x); 


außerdem ist infolge (11, 6) 
(25,5) Ana a)l Sf Tnla, ) O(WFL, z)|dt| = Aula, 2). 


Ist speziell 
(25, 6) IPe)| = pli)E, 
so wird infolge (24, 6) 


zT 





[pw ar) 
(25, 7) IY(a,z)| si + (e —1)E. 
Diese Ungleichung, zusammen mit 
Foo 1 
(25, 8) |Y(a, x) EsSe :E 


wurde von H. Fuhr [1, 2] angegeben. 
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Eine andere, von demselben Verfasser [1] herrührende Ungleichung ist eine Kon- 
sequenz der Formel (12,13). Nimmt man an, daß 


(25, 9) IP) < plt) E, |R()| <rü)E, 
also 

(25, 10) IE) = IPlt) + Ri] S (pt) + r(t)) E 
ist, hat man infolge (25, 8) 


t z 
| z [pe lar| + Sin) +rn)larl 
| |Y(a, x) — Z(a, z)| s fe ! .r(t) |dt|- E 


a 


(25, 11) 


fmerial freie 
\ = e® : (e? — 1) E 

Diese Ungleichungen können bedeutend verallgemeinert werden. Es sei P(x) in 
eine Relativmatrix 





mit quadratischen Diagonalelementen zerlegt, und jedes Element befriedige eine Un- 
gleichung von der Form 


|Pı|l Sgu Ew 


wobei g;, positive Skalare sind und 


Ei u Fi. Eım 
(25, 12) E -(: aus ) 
Emm 
also 
(25, 13) |Pals @uEa)- 


Wir nennen dann G = (g,,) eine Majorante der Relativmatrix P = (P;,). Unter dieser 
Voraussetzung wird 


(25, 14) W(a, x) = fa + Gt) | dt |) 


eine Majorante für die Relativmatrix | Y(a, x2)| —1 sein, d.h. 
) |Y;x(a, x)| s W «(a, x) Ex, i=k R 
| \Yı(a, x)) . 1 + (W;;(a, x) — 1)E; . 

Der Beweis hierfür beruht auf der Ungleichung 

E; E; < Ex 

(vgl. (4, 14)), die mit sich bringt, daß das w-te Glied in der Peanoschen Entwicklung 
für W(a, x) eine Majorante für das entsprechende Glied in der Entwicklung für Y(a, x) 
wird, falls dieses als Relativmatrix geschrieben wird. 

Mit Hilfe der eben bewiesenen Ungleichung kann auch (25, 11) verallgemeinert 
werden. Sind G,(x) und G,(x) Majoranten für P(x) bzw. R(x), und setzt man 


(25, 15) 


z 


[ m 
Wa, 2)= | (1 +6) dt) 
(25, 16) { 


Wıla, 2)= | (+ GC) + GW) de), 
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dann wird 
(25, 17) fwila, t) Gz(t) W;tt, x) | dt | = W;(a, x) —Wı(a, x) 


eine Majorante für Y(a,x)—Z(a, x) sein. Führt man eine Q-Transformation in W, 
aus, so bekommt man die direkte Erweiterung von (25, 11): die Majorante kann folgender- 
maßen geschrieben werden 


(25,18)  (J(1+ Wila,2)Gz) Wi, a) | de ))— 1) Wıla, 2). 


Bei Anwendung passender Q-Transformationen kann man oft nützliche Un- 
gleichungen erhalten. Hier soll nur eine einzige erwähnt werden, die wir gleich benutzen 
wollen. Wir nehmen eine Q-Transformation in Y(a, x) vor, mit 


zT 


Mm Fun a 
T= |U+Pd)=e 5 


wobei 


und bekommen durch (19,2) und (25, 4) }) 


B7 


| Ya, 2)|< | d + 17a, 2) IP) — Palı)| - TC, a)| dt) -|Tla, 2). 


a 


Da 
freue 
|T(a, 2)| = e* 
ist, erhalten wir dadurch, daß wir zurücktransformieren 
2 _ 
(25, 19) Ya, 2) s | + Pa), 
wo 
2 Rp, | Piel E Paul 
(25, 20) Ep 
|PaıbIPel Rdn 
ist ?). 


Wir berücksichtigen hier noch das unbestimmte Produktintegral. Da die Ab- 
leitung der Ungleichungen (25, 3, 4,5) und (25, 15) sich nur auf die Peanosche Ent- 
wicklung aufbaut, gelten die entsprechenden Ungleichungen auch für unbestimmte 


Produktintegrale, wenn nur Y(x) und W(x) ähnlich wie Y(x) gebildet werden. Be- 
sonders sei hervorgehoben, daß | Y (2)| s Y 4(2) und daß W;(x) bei geschickter Wahl 
von A eine Majorante für Y4(x) sein wird, wenn G(x) eine Majorante für P(x) ist. 


!) dt wird hier der Einfachheit halber als reell und positiv angenommen. 
*) Die Ungleichung kann auch unmittelbar aus der Definition des Produktintegrals abgeleitet werden: 


m—1 


m—1 FR 
IY(a, x)| = lim I |1-+ P(a,) öu| = lim IId+ P(a,) öu + o(öu)). 








us€s Pr as. 
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$ 26. Ein Grenzwertsatz. 


Zur Anwendung der gefundenen Ungleichungen wollen wir zuerst einen Grenz- 
wertsatz von Fuhr [1] beweisen: 

I. Wenn (P,(t)) eine Folge integrierbarer Matrixfunktionen ıst, die so gegen die eben- 
falls integrierbare Matrixfunktion P(t) konvergiert, daß 


z 


(26, 1) 1 P9— Pw|-Idt| 0, 


a 


dann wird 
N n 
(26,2) Y,la, x) = | (1 -- P,(t) dt) + (1 + Plt)dt) = Y(a,x) für v—w. 


Der Satz ergibt sich aus der Ungleichung 


[mc 1a! fear 
| Y,(a,2) — Y(a, 2)| s e* . (e® —A)E, 
wo 
IPW)ispt)&E, 1 PA)— P)| SA) E, 
die eine unmittelbare Folge von (25, 11) ist. 

Mit Hilfe dieses Satzes kann man eine eigentümliche Darstellung des Produkt- 
integrales, die von Schlesinger [3] herrührt, ableiten. Wird P(x) von stückweise kon- 
stanten Funktionen P/(t) approximiert, so kann Y(a, x) in eine gewisse Anzahl von 
Produktintegralen zerlegt werden, die infolge (24, 5) durch Quadratur bestimmt werden 
können. Für » > © ergibt sich der Ausdruck 

(26, 3) Y(a, x) = lim II ee, 


worin die Bezeichnungsweise analog der in $9 angewandten ist. 


$ 27. Gleichungen mit asymptotisch konstantem Koeffizienten. 


Seitdem die Poincaresche Abhandlung [1] erschienen ist, sind Gleichungen von 
der Form 

(27,1) Y'(z) = Y(z) P(x), 
wobei P(x) eine für 22 a> 0 stetige Funktion ist, die für x > © einem Grenzwert A 
zustrebt, oft untersucht worden. Da ich beabsichtige, in einer späteren Arbeit die hier- 
durch entstandenen Fragen in größerem Zusammenhange zu erörtern, soll hier nur ge- 
zeigt werden, daß man mittels der Ungleichungen einige wichtige Sätze von Perron be- 
weisen kann. 

Wenn P = A”' AA die Weierstrasssche Normalform von A ist, ersetzen wir Y durch 
YA; die Gleichung (27,1) wird dadurch auf die Form 

(27, 2) Y'(z) = Y(z)(P-+ H(x)), H(x2)-0 für 2 
gebracht. 

Wir beweisen zuerst eine Ungleichung, die für jede Lösung dieser Gleichung gilt. 
Es ist 


(27, 3) IY(z)| s |Y(a)| - |Y(a, &)|, 
und weiter nach (25, 19) 


|Y(a, 2)| S | (1 + MP + IH@))) di). 
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Es bezeichne nun o’ die größten der reellen Teile der Eigenwerte und n(x) } 0 sei so ge- 
wählt, daß 


(27, 4) |H(x2)| sn(z)E. 
Indem weiter 
00...00 
(27,5) o-(19. 2! 
00. 10 


folgt hieraus 
na ae 
(27,6) |Y(a, 2) < [1 + (e + © + nl) E)dı) = e"""7 la, 2), 


mit 
z 


(27,7) V(a, x) = fa + (© + nit) E) dt). 


a 


Zur Abschätzung dieses Produktintegrals nehmen wir eine Transformation vor: 
_ nn 
Vla, 2) = e® / (+ eEe* (1) dt) 


und setzen die Peanosche Reihe für das neue Produktintegral ein. Dadurch bekommen wir 
(27, 8) ”(a, 2) = 2V,(a, 2), 


wobei 


(u) 


(27,9) V,(a,2)= fe WI EeMWE,,. EM Walt)... nlt)di... dir. 


a 


Nun läßt sich offenbar zu jedem positiven e eine Konstante y derart bestimmen, daß 
(27, 10) e®"<ye”E für >00. 

Folglich gilt die Ungleichung 

(yf nt) di) 


<E, 
u! 





\ *u) 
V.(a, 2) <a)... nltu)dt,... di E= ye 
und hieraus ergibt sich 


z 
z—a)+ [me)dt 


(27, 11) V(a,x) < ye “a  .E. 
Mit beliebig kleinem ö ist also 
(27, 12) |Y(a, x)| < ce +9eo E, 


wobei c von ö, nicht aber von x und a abhängt !!). 

Durch Anwendung der Formel (16, 7) sieht man, daß 

(27, 13) |Y(z, a)| < ce + 5, 
wo 0” der kleinste unter den reellen Teilen der Eigenwerte von P ist. Hieraus kann man 
eine untere Grenze für gewisse Elemente in Y(a, x) ableiten. Da 


1 <|Y(a, 2)| - |Y(z, a)| < ce +9. |Y(a, x)| E 


ı) Perron [2], S. 28. 
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und 
1 <|Y(z, a)| -|Y(a, z)| S ce” +9. E|Yl(a,)l, 
hat man nämlich sowohl 


’ n "—dXz—a 
(27, 14) ei Y.(a, z)| 2 2 7 
als auch 
(27, 15) | &]y,(a, 2)| > a ee 


Es gibt also in jeder Zeile und in jeder Spalte von Yf(a,x) mindestens ein Ele- 


ment > — ee Nee) A), 
c 

Für spezielle Lösungen lassen sich viel bessere Ungleichungen ableiten. Wir setzen 
zuerst ‘voraus, daß die reellen Teile r,,...r„ der Eigenwerte von P ungleich sind und 
zwar, daß 

(27, 16) H>R> ++ > 
Wir wollen dann die Lösung 

(27,17) Y(z) = a" Vı(z) 


betrachten, bei der 
F- z 
(27,18) Vı(z) = / 1 + ee Puı)e"9a)=e” / (1 -- eH(t) e”" dt) e* 
(L) (L) 


(vgl. $ 15) wobei 


aar-.:-ada 
(27, 19) | 
oo » a 
Wird 
5 r, 0 ...0 
(27, 20) -w-( Br 0 | 
DI - +. Mm 


gesetzt, dann gilt nach der Schlußbemerkung in $ 25 die Ungleichung 


z 
rn 


(27,21) |Vı@a)! s J + et ER" zu) dt) = Vıle). 
(2) 
Die Glieder der Peanoschen Reihe für V,(x) sind durch die Rekursionsformel 


(27,22) Turtle) = Se" Tue) E” nd) Jar], Tale) = 1 
(L) 
bestimmt (vgl. (19, 9)). Um eine brauchbare Ungleichung für V;(x) hervorzubringen, 
suchen wir die skalare Funktionenfolge ( p,(2)) derart zu bestimmen, daß 


(27, 23) T(@)<p,()E n>21, la). 


Es muß dann 


Ta) Sf e "Er" ot) mt) ld. 
(L) 


Die Elemente der Matrix rechts sind teils 





1) Vgl. Perron [1], S. 26. 
Journal für Mathematik. Bd. 171. Heft 2. 13 
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z 


h$  -nt-2) < 1 u 
fe pt) nl) dt S si fr, p,(t) n(t) dt für ı Sk, 


a 


teils 
an 


BEER N en 
ran as, [em n,0m@d tür i>k, 


z 


wobei ö =Min|rr —r;|. Es Fr also 
(28, 24) Pur? -/ 9, O)mtı) di + j Ep, ml) de 
zu wählen. Zur Untersuchung dieser Folge führen wir die Funktion 


(27, 35) Pla; A) = &p,(a)i* 
ein. Sie befriedigt die Integralgleichung 


(27,56) pa; )=1+Af PU; A) nt)d+ fe" "pl; And) dt, 


woraus man leicht die Differentialgleichung 


(27, 27) 9” (2; A) = 6p'(x; A) — Aön(x)p(x; A) 7 
ableitet. Für | 
(27, 28) o(z; A) = (log p(x; A))' 
gilt dann die Riccatische Differentialgleichung 
(27, 29) o’(2; 4) = Öolz; A) — 0?(x2; A) — Aön(e). 
Die Glieder der Potenzreihe 
(27, 30) o(2; 1) = o,(2)2 


sind bestimmt durch die Rekursionsformeln 


0,2) = 6 fer 64_,(t) o,(t) dt 
(27, 31) { ä 





=ödfje* E Ou—(X +) o,(z + t) dt 
0 gr 
und 
(27, 31a) o,(2) =d fe Pnl)d = öde tn +1) di. 
z 0 


Durch Induktion erkennt man, daß jedes o,(x) monoton abnimmt, und hieraus er- 
geben sich die Ungleichungen 


(27, 32) o,(2) <o,(2), v1, 
wobei 

| 4 

(27,33) la) = nla) und 5,2) =8,5,-,)öe), 22. 

Die Funktion 

















(27, 34) (x; A) = 2a,(z) A 
befriedigt die quadratische Gleichung ; 
(27, 35) o(z; A) = (o(z; 4))’ + An(e), 4 
so daß. } 
2 1 ' i 
(27, 36) o(2;,))= DRL — 1 — 4/n(2)). | 








ie RR 


EEE 
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Da o(x; A) für | A| <s eine reguläre analytische Funktion von A ist, wird 


auch o(z, A) und p(z; 4) regulär in diesem Kreis sein. Damit haben wir die Existenz 
von V;(x) unter der Voraussetzung | 


1 
(27, 37) n(x) < zT 
bewiesen, und weiter ergibt sich die Ungleichung 
zw 
(27, 38) |Vılz)| se -E, oft) = oft;1). 
— Es sei übrigens bemerkt, daß 
(27, 39) (2) = (2) + O(n?(z)) 
und somit 
(27, 39a) fett) di = [ntt) dt - (1 + o(1)) = o(e). 
Ganz analog beweist man für 
(h) 
(27, 40) Vz!(2) = / (1 — ED ul) ee-® an) 
(vgl. (15, 12)) die Ungleichung 
Fa ae 
(27, 41) IVve@alse -E, 
woraus man schließt, daß es sowohl in jeder Zeile als auch in jeder Spalte von V;(r) 
fra 
mindestens eın Element > gibt. 


Der allgemeine Fall, bei dem keine Voraussetzungen über die Eigenwerte von P 
gemacht werden, läßt sich mit fast derselben Methode erledigen. 


(27, 42) -( 
+ + 


wird so geordnet und eingeteilt, daß die Eigenwerte von P; denselben reellen Teil 
r; haben und weiter 


(27, 43) n>Nn>*''>Tn 
Wir betrachten dann die Lösung (27, 17), wo diesmal 
(27, 44) L= (Lu), 


und Z,;; eine quadratische Matrix von derselben Ordnung wıe P; ist, und die Elemente 
von L;« gleich a oder & sind, je nachdem ı Ss oder > k ist. 


Wird 
P‚ r, en; RM 2 
(27, 45) -( En und RP=P--© 
Fe 
gesetzt, so genügt das Polynom e9=Ee®% für x > 0 und beliebiges e > O0 der Ungleichung 
(27, 46) e9% Ee9% < ye=E, 


wenn y = y(e) passend gewählt ist. Es ist also 


zT 
MM 


(27,47)  |Vıle)! < Vılz) = / (1 -- yert-» Kee-9 eelz-t „(t) Idt|). 


(L) 
13* 
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In den vorhergehenden Erörterungen ersetzen wir also einfach n(t) durch ye’"!n(t). Die 
(27,26) entsprechende Integralgleichung lautet dann 


(27,48) plz; A)=1+ Ay fe@ dal; A)nlt)di+ Ayfert Mm gl; Anl) dt. 
Die Funktion y(xz; A) = e”"o(x; A) befriedigt die Gleichung 
(27, 49) y(z; 4) = e7= + Ay fylt; Aynlt) di + Ay fer ea ylt; A)n(e) dt 


und hat wegen e-"”" < 1 die oben gefundene Lösung von (27, 26) als Majorante, wenn / 
durch Ay und ö durch ö—.2e ersetzt wird. 
Hieraus folgt 


(27, 50) | Vı(z) | <etete E. 
Dieses Ergebnis läßt sich aber sofort zu 
(27, 51) |Vılz)| < et® E 
verschärfen. 
Hierzu bemerken wir erstens, daß — mit naheliegenden Bezeichnungen — 
(27, 52) Vralr) = e?*C(a, b) e®* Vinyl), y 
wobei die Konstante | 
„0 
(27, 53) C(a, b) = / (1 + etH(t)e de) 
L(a) 
sich als eine relativ halbreduzierte Matrix erweist. Sie ist nämlich als der Grenzwert 
L(b, e) | 
(27, 54) lim (1 -- efH(t) e="' dt) 


e>n 


L(a, e) 
aufzufassen, indem 
aı Sk 
Lia, ) = (Lu(a, ), Lula,) = 10 3% 
Wird nun H(t) ähnlich wie P eingeteilt, 
H(t) = (Haft), 
so läßt sich der Produktintegrand in (27, 54) in der Form 
1 + (ef Hult) e’# dt) 
schreiben [ vgl. (15,8) und (15, 7)], und hier muß offenbar füri>k dti= 0 gesetzt werden. 
Bezeichnet H(t) diejenige Matrix, in die sich H(£) verwandelt, wenn jedes H;.(t) mit > k 
durch O ersetzt wird, so ist also das Produktintegral in (27,54) unabhängig von c und gleich 


RT Fans VRETE Fe 





h 
Mm — 
(27, 55) C(a,b) = [ (1 + et Alt) e-® dt) = (Cu(a, b)) 


mit C;x(a,b) =0Ofüri> k. Der Faktor von V;«u,(x) in (27, 52) ist somit 
(e""* Cala, d) e*) = O(2), 


wobei x eine passend gewählte ganze Zahl ist. 
Weiter bemerken wir, daß die Herleitung der Ungleichung (27, 50) nach (27, 37) 


1 ) 
yn(a) < % voraussetzt. Nachdem e gewählt ist, bestimmen wir also ein 5 derart, daß 


yn(b) < 7 ist; aus (27, 52) und (27, 50) ergibt sich dann 
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(27, 56) Vxcay(2) = Ole“) 
für beliebig kleines e. Bezeichnet v;.(x) die Elemente in V;(x) = Vı«,(X), so hat man also 


(27, 57) Tim 8 | in <o0 


und aus dem entsprechenden Ergebnis für V7'(x) schließt man weiter, daß es sowohl in 
jeder Zeile als auch in jeder Spalte mindestens ein Element gibt, für das 


(27, 58) tion _0. 


In bezug auf die Anwendungen ist es nützlich zu bemerken, daß es für unsere Über- 
legungen durchaus nicht wesentlich ist, daß P gerade eine Weierstraßsche Normalform 
ist, es genügt, wenn es eine halbreduzierte Matrix ist. Die einzige hinzukommende 
Schwierigkeit rührt daher, daß © nun nicht mehr permutabel mit P zu sein braucht. Man 
kann sie aber dadurch umgehen, daß man in der Definition (27, 18) von V,(z) die Ex- 
ponentialfunktion e®* durch eine passend gewählte Lösung der Gleichung 


(27, 59) &’ —=&P 
ersetzt. Man wähle z.B. 


x 
MN 


(27, 60) Plz) = er. / (1 + efrtt-2) Hefriz- dt) = efr Y,(x), 
(h) 

wobei P, diejenige Diagonalmatrix bedeutet, die dieselbe Diagonale wie P besitzt, und 
L durch (27, 44) definiert ist. Da die Peanoschen Reihen für Yj'(z) nur eine endliche 
Anzahl von Gliedern enthalten — man vergleiche den Schluß von $ 24 — weist man 
leicht nach, daß diese Funktionen O(z”-!) sind, und das genügt für unsere Überlegungen. 

Ist nun eine skalare Gleichung r-ter Ordnung mit asymptotisch konstantem Ko- 
effizienten 


(27,61) ym + (a + la) rd + + (an + nlz))y = 0 
gegeben, so ersetzt man sie durch das Gleichungssystem erster Ordnung 





(27, 62) ae en 
Ym 91 t U 
su trmht tm 
wobei die n, lineare Verbindungen von den e,sind und 9,,... 0, die Wurzeln der determi- 
nierenden Gleichung 
(27, 63) Tu Sk 


bezeichnen. Da die Differentialquotienten von y leicht angebbare lineare Verbindungen 
der y, sind, ergibt sich sofort aus (27, 57) und (27,58) der Perronsche Fundamentalsatz 
[1], S. 46: 


I. Es existiert ein Fundamentalsystem y,,-..y, zu (27,61), derart, daß für jedes i 
erstens 

(27, 64) en u are. 
und zweitens rg 

(27, 65) lim nu base — Ro, 


für mindestens ein k<n—1. 








102 Rasch, Zur Theorie und Anwendung des Produktintegrals. 


Kapitel III. 
Gleichungen mit analytischen Koeffizienten. 

$ 28. Mehrdeutigkeit der Lösungen in einem Bereich. Windungsfunktionen. 

Im folgenden soll die Funktion P(x) nicht nur auf einer Kurve £ gegeben sein, 
sondern in einem Bereich Q, in dem sie eindeutig und regulär analytisch ist. 

Ist A einfach zusammenhängend, so sind die Glieder in der Peanoschen Reihe für 
Y(a, x) regulär und nur von den Integrationsgrenzen abhängig und nicht von der spe- 
ziellen Kurve in X, die sie verbindet. Hieraus folgt der Satz (Volterra [2], S. 22]).: 

I. Wenn P(x) in einem einfach zusammenhängenden Bereich regulär ist, dann wird 
Y(a, x) es auch sein, und ist & eine beliebige geschlossene Kurve in seinem Innern, so wird 


(28, 1) fa + Plt)dı) =1. 


Wenn X nicht einfach zusammenhängend ist, dann wird Y(a, x) im allgemeinen 
mehrdeutig sein. 

Ist x’ ein Punkt über x auf der Riemannschen Fläche über W, so kann man den 
Integrationsweg von a nach x’ infolge des eben bewiesenen Satzes durch den Punkt a’ 
legen, der in derselben Weise über a liegt, wie x’ über x. Hieraus folgt 


Y(a,x©) = Y(a,a’)Y(a’, x’). 


Setzt man 

(28, 2) U(z) = Y(z, x) 
und beachtet, daß aus der Eindeutigkeit von (x) 

(28, 3) Y(a’, 2’) = Yla,:) 
folgt, so hat man die Beziehung 

(28, 4) Y(a, x’) = U(a)Y(a, x), 


welche die Mehrdeutigkeit von Y{a, x) charakterisiert. 


Betrachtet man eine beliebige Fundamentallösung Y(x), so ıst Y(x’) = Y(x) 
offenbar auch eine Fundamentallösung. Die Konstante C in der daraus folgenden Be- 
ziehung 


(28, 5) Y(x«)=CY(x) 
kann dadurch bestimmt werden, daß man Y(x) durch Y(a, x) ausdrückt; man be- 
kommt dann 

(28, 6) C = Y(a)U(a)Y”(a). 

In einem gegebenen Bereich lassen sich im allgemeinen unendlich viele Funktionen 
von der Form U(x), die wir Windungsfunktionen nennen wollen, definieren. Diese Funk- 
tionen bilden offenbar eine Gruppe, die aus einfachen Windungsfunktionen zusammen- 
gesetzt werden kann. Jede Kurve von endlicher Länge auf der Riemannschen Fläche 
über W, die x mit einem Punkte über x verbindet, kann man nämlich deformieren und 
dann in eine Reihe von Kurven auflösen, deren Projektionen auf die x-Ebene jede eine 
Windung um eine Randkurve in X beschreiben. Nennt man die Windungsfunktion, die 
die v-te Randkurve einmal in positiver Richtung umschlingt, U,(x), so wird die betrachtete 
Windungsfunktion in Faktoren U#'(x) aufgelöst. Speziell sei bemerkt, daß für einen 
doppelt zusammenhängenden Bereich alle Windungsfunktionen Potenzen von einer ein- 
zigen sind, da die Windungsfunktion, die p Windungen um die innere Randkurve des 
Bereichs entspricht, die p-te Potenz derjenigen ist, die einer Windung entspricht. 

Die angeführte Gruppe von Windungsfunktionen bildet einen Teil eines Funktionen- 
körpers. Die Windungsfunktionen befriedigen nämlich infolge (10, 6) die Differential- 
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gleichung 

(28, 7) U’(x) = U(z)P(x) — P(a)U(e), 
die infolge (20, 11) die Lösungen 

(28, 8) U(x) = Y’(a, x)C.Y(a, x) 


hat, wo C alle Konstanten durchläuft, und hieraus folgt, daß Summe und Produkt zweier 
Lösungen wieder eine Lösung ist, sowie, daß die reziproke Matrix zu einer invertierbaren 
Lösung auch eine Lösung ist; schließlich wird die Gleichung auch von der Einheitsmatrix 
befriedigt. 

Hieraus folgt unter anderem, daß auch (o — U(xz))”" eine Lösung sein wird, wenn 
U(x) eine beliebige Lösung und o ein skalarer Parameter ist. Weiterhin sieht man, indem 
man (o — U (x))”" in Partialbrüche zerlegt, 


(28, 9) (e — U(2))” “ar ep: 


daß die Funktionen U ,,(x) die Gleichung a weil die Eigenwerte von 
U(x) unabhängig von x sind. 





$ 29. Die Darstellung der Lösungen in einem doppelt zusammenhängenden Bereich. 


Wenn X doppelt zusammenhängend ist, wird die Mehrdeutigkeit einer Fundamental- 
lösung vollständig charakterisiert durch eine Beziehung 


(29, 1) Y(a)=CY(e), 
wo x’ über x auf der Riemannschen Fläche über W liegt. Setzt man 
1 


wobei man den Logarithmus beliebig wählen kann, dann wird x mehrdeutig von derselben 
Art sein, vorausgesetzt daß x = 0 innerhalb der inneren Randkurve von W liegt. Daraus 
folgt, daß x=4Y(x) eindeutig in Yist. Hat A die Normalform P, so gibt es eine Funda- 
mentallösung, die man 

(29, 3) Y(z) = z’V(z) 
schreiben kann, wo V(x) eindeutig in W ist. 

P ist im Vorhergehenden durch den Wert einer Windungsfunktion in einem ein- 
zelnen Punkt bestimmt. Um weiterhin V(z) eine solche Form zu geben, daß deren Ein- 
deutigkeit in A sofort zu erkennen ist, machen wir Gebrauch von einem Kunstgriff, der 


im Prinzip von Frobenius [1] herrührt. 
Da x’ = x e?”i ist, zeigt man durch direkte Ausrechnung, daß 


(29, 4) Y(z; 0) = ft" Yıt, x’) dt - (e — U(x))" 


eine Lösung der inhomogenen Gleichung 

(29, 5) Y'(2) = Y(z) Pix) + a7', 
ist, wobei o einen komplexen Parameter bezeichnet. Bei der Ausrechnung kann man be- 
nutzen, daß (e”"* — U(x))”" die Differentialgleichungen (28, 7) der Windungsfunktionen 


befriedigt. 
Um Lösungen der homogenen Gleichung zu finden, sei bemerkt, daß die Eigenwerte 


0,» .».». ® Om von 


(29, 6) U(z) = — ; log U(x) = 2 7 "(a x) log U(a) - Y(a, &) 


. unabhängig von x sind, be daß Y(z; o) als Funktion von o meromorph ist, mit den Polen 
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je1l..M 

„=d, +1... 

daß das Residuum von Y(x; 0) in einem beliebigen Pol eine Lösung von 
(29, 7) Y'(z) = Y(z) P(x) 


in den Punkten 0, + », oder doch in einigen davon. Hieraus folgt, 


ist. 

Dabei machen sich nun zwei bedeutungsvolle Verhältnisse geltend: Einerseits kann 
man mit einer endlichen Anzahl dieser Partikularlösungen eine Fundamentallösung auf- 
bauen, und andererseits ıst jede einzelne Lösung ein Produkt von einem x°' und einem Poly- 
nom in log x, dessen Koeffizienten in A eindeutige Funktionen sind, so daß man eine Lösung 
von der Form (29, 3) konstruieren kann. 

Um das erste zu erkennen, berechnen wir 

m m 

(29, 8) Y,(z) a Res. Y(x; oe) Y,(z). 

g=0;trV vn 
Hierzu sei einmal bemerkt, daß man (29,4) auch 


(29, Aa) Y(z; o)=[ Yit, x) PR TREIR 
schreiben kann, und dann, daß man bei passender Wahl von U(z) — nämlich so, daß 
0,1»: O,„, keine ganzzahligen Differenzen #0 haben — für die Bildung der Residuen- 


summe Y,(z), den Faktor des Integrals durch —i (o— v»—U(x))”" ersetzen kann. Nach 


(7,1) hat man also 
Br ä Wz)+r—1 
(29, 9) Y‚(z) m) Y(ı,2)i dt. 


Nun weiß man ja im voraus, daß die Gleichung (29, 7) eine Lösung von der Form (29, 3) 
hat, und daraus schließt man leicht, daß mit passend gewähltem A 


(29, 6a) Ulz) = Y’'(z)PY(x) 
ist. Somit haben wir 


’ 


1 f —1 P+v—1 1 [ —l v—1 
a en = —. — 
(29,10) Y,(x) mi) Y(t)t dt- Y(x) Im V"()t dt-Y(x) = B,Y(x), 


wo B, der Koeffizient zu a” in der formell gebildeten Laurentreihe in X für V"(a) ist. 
Diese Reihe ist nach Mittag-Leffler !) wenigstens in Y summierbar, indem z.B. 


en . = B,a’ 
(29, 11) V” (a) = lim ri BES... EEE 
= re) 
ist. Wie man leicht einsieht, kann man diesen Ausdruck durch 
ih . 2 B a” 
(29, 11a) V” (a) = lim nn, 
k>xT Pe + ı) 


ersetzen, und da V "'(a) invertierbar ist, muß der Ausdruck unter dem lim-Zeichen es auch 
sein, wenn k hinreichend groß gewählt ist (vgl. $5, Schluß). Aber hieraus folgt, daß, mit 


solchem X, 
+k! 


% Y,(z) Fam ) 





1) Atti del IV Congresso Mat. Roma 1908, S. 82. 
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eine Fundamentallösung ist, was mit sich bringt, daß die Sammlung von den Matrizen 
Y.(@),Jj=14...m,v=0, +1,...+Kln linear unabhängige Vektorlösungen enthalten. 


Durch Anwendung der dritten seen Form (8, 10) von t"” bekommt man 


U;.(z) 
ja! 


h) 


(29,12) Y„(e)= aut”. ef Y(zt, x) ii’! (log t -- log x)" dt: 


indem 
(29, 13) e-W- FI" mr 


Da nun sowohl U,,(x) als auch 
1 


Mm 
(29, 14) Ylzt, x) = / (1 + zP(zt) dr) 
t 
eindeutig in X sind, sieht man, daß der Faktor zu x® ein Polynom in log x ist, wobei die 


Koeffizienten eindeutig in W sind. 

Wir betrachten nun die Vektorlösungen in Y,,(2), v»=0, +41,...+ kl. Wir wählen 
darunter eine aus, deren Faktor von x°’ vom höchsten Grade in log x ist. Den Grad nennen 
wir y,— 1, und die Lösung schreiben wir in der Form 

Yıml 


log x)” 
Ya, (2) = 1,1%) u El 


Man kann nun sofort y,, — 1 andere Lösungen angeben, da die Funktionen 


u—1 2 


29,15) DR hen 


linear ausgedrückt werden können durch Y,,(®) und deren analytische Fortsetzung 
zu den Punkten x: e=?,p= +1, +2,... und folglich unsere Differentialgleichung be- 
friedigen. 

Die Sammlung von Vektorlösungen, die zu x°’ gehören und von der wir ausgingen, 
wird nun reduziert, indem wir y,, passend gewählte Lösungen durch das System (29, 15) 
ersetzen und die übrigen mittels dieses reduzieren, so daß die Polynome in log x von so 
niedrigem Grad wie möglich werden. Ausgehend von einer Lösung, die danach log x in so 


hoher Potenz wie möglich enthält, bilden wir nun ein Lösungssystem analog zu (29, 15): 
‚u 


4 
(29, 16) Y.(%) - In, u j* (og: d „ u=h... Ya‘ 


Mit Hilfe der gefundenen zwei TEEN - 15) und (29, 16) wird die Reduktion 
fortgesetzt, indem noch ein Lösungssystem konstruiert wird, und das wird solange fort- 


gesetzt, bis die Reduktion vollendet ist. 
Die gefundenen Vektorlösungen lassen sich zu einer rechteckigen Matrixlösung 
von einfacher Struktur zusammensetzen. Führt man die Bezeichnungen 


B,,(@) 00...00 
(29, 17) 1.1 =| ZERE ) -[10.. EZ )ı- (0) für y„=1) 


ein, so kann man schreiben 
(29, 18) Y.(x) = 


Journal für Mathematik. Bd. 171. Heft 2. 14 
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und setzt man weiterhin 
61 +0 V;ı(®) 
(29, 19) P,=g,+6,, = x Be V,(x) = eb, 7 
%; j 


so hat man 


1 (x) 
(29, 20) Y (x) -( ie } == ZH Vulz). 
Y„,(&) 

Es ist klar, daß sämtliche zum Exponenten go, gehörigen Lösungen sich linear aus 
den Vektorlösungen in Y,;(x) bilden lassen. Wir müssen aber noch nachweisen, daß diese 
Lösungen linear unabhängig sind. Nach (29, 8) und (29, 10) ist jede Vektorlösung in 
Y,;(x2) von der Form 


(29, 21) D,,(a) = a7 B,,-zV,(a) = zo zer g 2)" 8,0 Ve), 


wobei V;(x) eindeutig in X und ; ähnlich wie ©; gebaut ist, und hieraus schließt man 
ohne Schwierigkeit, daß der geschilderte Prozeß wirklich zu linear unabhängigen Lösungen 
führt. 

Die m Partikularlösungen Y;(x) kann man weiter zusammensetzen zu einer Fun- 


damentallösung 


Yı(z) 
(29, 22) Y(x) -( Kr ) 
m(%) 
und diese ist genau von der Form (29, 3), wenn wir 
WE Vı(2) 
(29, 23) -| EEE ER )) V(x) -( e 
0... P Vn(z) 


setzen, und damit haben wir unser Ziel erreicht. 
Im übrigen kann man durch effektivere Ausnutzung der Normalform von U(x) 


eine bedeutende prinzipielle Reduktion des Darstellungsproblems erreichen. 

Da 
| U(z) = Y'(x) e" Y(z) 

— VU(x) &" V(x) 

ist, kann man, wenn U(x) gegeben ist, durch Lösung von linearen Gleichungen eine in 
A eindeutige Funktion 7(x) bestimmen, für die 

(29, 24a) U(xz) = T'(z) e”* T(x) 
ist, und transformiert man das Produktintegral Y(a, x) damit, so bekommt man 


(29, 24) 


(29, 25) Y(a, x) = 1a) [di + MP] T)dt) T(x). 


Da U(x) die Differentialgleichung (28, 7) befriedigt, erkennt man durch Anwendung 
der Transformationsformel (20,5), daß @ = Q(P|T) permutabel mit e" ist und folg- 
lich auch mit P selbst (vgl. $7, I). 

Wenn nun P,,... P„ in (29,419) von der Form 


00::-.00 
(29, 26) mat er Bi: 
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sind, dann liegt die Situation vor, die in Volterras Satz über permutable Matrizen ($ 3, I) 
beschrieben ist, und daraus folgt, daß Q(t,) und Q(t,) permutabel für beliebige t, und t, 
in X sind. Aber dieses führt mit sich, daß die Produktintegration über 1 + Q(t) dt durch 
Quadratur (vgl. $ 24, 2)) ausgeführt werden kann, also 


f Qıyae 
(29, 27) Y(a, x) = T"(a) e“ : T(z). 
Das Residuum 
1 P7 


ist augenscheinlich auch mit P permutabel, also auch mit @(t) und weiterhin mit 


(eo) dt, so daß die Beziehung (8, 3) angewandt werden kann. Damit haben wir 


(29, 29) Yla, 2) = T*(a)-(7) - vie), 
wobei 


Fan-A 
(29, 30) V(x) = e - T(x) 
eindeutig in W ist. 
Aber selbst wenn dieser günstige Fall nicht eintritt, kann man bedeutende Vor- 
teile durch die Ausführung der Transformation (29, 25) erlangen. Transformiert man 
nämlich weiterhin mit x’ und erinnert sich, daß xP permutabel mit Q(t) ist, so bekommt 


man 


(29, 31) Te + Odt) = fu +01 )dı) a — (£) | fu + (e——) di) 
und hier ist 
(29, 32) Zia, 2) - {1 + (a) di) 


a 


eindeutig, da die Windungsfunktion Z(zx, x’) offenbar = 1 ist. Um eine Darstellung von 
Z(a, x) zu erhalten, aus der die Eindeutigkeit unmittelbar hervorgeht, kann man so vor- 
gehen wie bei der Konstruktion von Y(x). Der Prozeß ist aber dieses Mal wesentlich 
einfacher, da alle Pole von Z(z; o) einfach sind und in den Punkten 0, +1, + 2,... 
liegen; wir haben also nur n linear unabhängige Vektoren aus den Matrizen 


g* 
u Ri. 1 [ 7 BE 
(29, 33) ZA) = x Zi. Z(at, 2) U di 


auszuwählen. 
Es tritt außerdem noch eine weitere Vereinfachung ein. Da Q permutabel mit P 
ist, muß es nämlich von der Form 


6, rg 
(29, 34) e-| Bern E 
ia 


sein, wobei Q; von gleicher Ordnung wie P; ist, mit dem es im übrigen auch permutabel 


ist. Z(a, x) ist also auch eine relative Diagonalmatrix, deren Diagonalelemente 
14* 





108 Rasch, Zur Theorie und Anwendung des Produktintegrals. 


n P, 
(29, 35) 2, 2)= (1 + (9) a) 
sind, und Z,(z) wird in entsprechender Weise in 
Z,ı(z) wie 0 
(29, 36) Z,(x) -( Bee 
0 ++ Zumlz) 
eingeteilt. 


8 30. Die Laurentreihe für Y(x). 


Bei der Berechnung von Y(z; o) und Y,(x) wird man wohl oft mit Vorteil die im 
allgemeinen rasch konvergierenden Peanoschen Reihen für Y(zt, x) und U(x) benutzen. 
Wenn es aber von Interesse ist, die Laurentreihe für V(x) aus der Laurentreihe 


(30, 1) P(x) = z A,ar-ı 


zu bilden !), so kann das auf folgende Weise geschehen: 
Das Produktintegral 


fi 
(30, 2) Ylz, iz) = / (1 + zPiat) dt), 


das für A = e?” in U(z) übergeht, ist eindeutig in X und hat also eine summierbare 
Laurentreihe 





(30, 3) Yiz, iz) = ZB.) 2* 
indem z.B. 
XS Bi) 
(30, 3a) Y(x, Az) = lim Zro]:| 74) 


gleichmäßig in einem beliebigen geschlossenen Bereich X im Inneren von X gilt. Setzt 
man nun 





A 
nn Pl) EN 
so wird infolge $26, I 
2 
(30, 5) Y.(z, Az) = / (1 + zP,(zt) dt) > Y(x, Az) 
1 


in. Da ?,(x) keine anderen Singularitäten als 0 und © hat, wird dasselbe für Y,„(x, Ax) 
gelten, so daß diese Funktion eine konvergente Laurentreihe 
+0 


(30, 6) Yılz, Ar) = 3 Bu(A; m) 2° 


hat. Die Koeffizienten hiervon werden dargestellt durch 
1 au 
(30, 7) Bulasm) = u: | Yale, a) 0" ar, 
g 


wobei & eine geschlossene Kurve um 0 ist. Wird sie in X gelegt, so erkennt man aus 
(30, 5), daß 
(30, 8) lım B.(}; n) = B,(A) 


n>0 


1) Dieses Problem wurde zum ersten Mal von H. v. Koch [1] gelöst, als eine schöne Anwendung seiner 
Theorie der unendlichen Determinanten. 
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ist, und damit sind die Koeffizienten in der summierbaren Reihe (30, 3) selbst dargestellt 
durch ein gewisses Summierungsverfahren. 

Analog hierzu werden die Laurentreihen, die wir im folgenden aufschreiben, summier- 
bar sein und ihre Koeffizienten werden meistens durch ein Summierungsverfahren analog 
dem für B,(A) beschriebenen bestimmt werden. Wenn wir schreiben, daß ein Koeffizient 
durch eine unendliche Reihe bestimmt ist, meinen wir deshalb nur, daß die Reihe summier- 
bar ist, in dem Sinne, daß der Koeffizient gleich ist dem Grenzwert für 7 — 0 der Reihe, 





die entsteht, wenn jedes A, durch ersetzt wird. 


A, 
n|x!+1) 
Für das w-te Glied in der Peanoschen Reihe für Y(z, Az) gilt offenbar 


+» 
(30, 9) T,(2; 4)=& B£(i)a‘, 


wobei die Koeffizienten durch 
(30, 10) BP(}) = Du. Anr A 


” ....H =“ 
+ u 


[2 


bestimmt sind und 
fi 
(nu) u 


(30, 14) 7 EFT © 7 5 2 un RHRE Se Re = 
1 


ist. Um diese Integrale zu berechnen, sei bemerkt, daß ein w-faches Integral als Funk- 
tion von 4 eine lineare Differentialgleichung w-ter Ordnung befriedigt, die sich leicht in 
eine Gleichung mit konstantem Koeffizienten transformieren läßt. Löst man diese nach 
einer bekannten Methode, die von Cauchy herrührt, so findet man, wie Günther [1], 
S. 302—303 gezeigt hat, folgenden einfachen Ausdruck 


Fr 


A 
(30,12) U; ru...) =— Den ee: 


Danach sind die Koeffizienten in (30, 3) durch die konvergente Reihe 








(30, 13) B,(i) = & BY(9) 
.- 
bestimmt. 
Wird x und A durch xt und t”’ ersetzt, so bekommt man die Laurentreihe 
+ @ 
(30, 14) Y(xt, x) =& CA), 
deren Koeffizienten durch 
(30, 15) a) = 3"), 
(30, 16) cd) = PP et; 1... 4) An... Any 


“tr tr, =y4 


’ E 
RR A Sa Ä E E 
bestimmt sind. 

Um nun z.B. Y,(x) zu berechnen — Y,;(z) und Y(x;o) kann man ähnlich be- 
handeln — sei bemerkt, daß die Integration in bezug auf tin (29,9) unter dem Residuums- 
zeichen in (30, 47) ausgeführt werden, kann auf Grund der gleichmäßigen Konvergenz 
und Summierbarkeit in bezug auf t von (30, 15) bzw. (30,16). Wir bekommen dann 
den folgenden Ausdruck 

+» © 


(30, 18) Y,(x) = & äh, ... A,,2(U(z) 7 9 a, '°* Xu) u U) 


© n=0 
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wobei 
2rilu—.) ___ 1 1 
30, 19 3 nn — —, 
( x ’ .) io 2nulu —L) SE+ 2%)... (+++) 
Um zu der gewünschten Darstellung von Y,(x) zu kommen, muß man hier z. B. 
die dritte kanonische Form für die auftretenden Funktionen von U(x) einführen, und da- 
zu braucht man die Laurentreihen für die Funktionen U;,(z). Um diese zu bestimmen, ist 


zu bemerken, daß sie durch die Funktionen U,;,(x) ausgedrückt werden können, die durch 





(30, 20) (wo — U(z))" = 22 u _ ET 


definiert sind, da die dritte kanonische Form für (o — U(x))”" = 2ri log” (wU”"(x)) 
offenbar 


a er 

(30, 21) (o — Ulx))" = zu o*(e “ 9,) 

u: 0w; 
ist. 

Da wir die Laurentreihe für jedes Glied in der Matrix ® — U(x) kennen, können 
wir die Laurentreihe für deren Unterdeterminanten (rn — 4)-ter Ordnung und damit 
für die komplementäre Matrix (v — U(x))”' - Det. (o— U(x)) bilden. Die Partial- 
bruchzerlegung von (vw — U(x))"" kann dann vorgenommen werden, wenn 

(30, 22) o(®) = Det. (wo — U(x)) 
in Faktoren zerlegt ist. 


p(w), das unabhängig von x ist, kann auf verschiedene Weise gebildet werden. 
Man kann U(x) in einem beliebigen Punkte in X berechnen und dann die Reduktions- 
determinante in diesem Punkt ausrechnen. Man kann auch von der Laurentreihe 
für U(x) ausgehen, da nämlich 9(w) das konstante Glied in der Laurentreihe für 
Det. (vo — U(zx)) ist. Eine dritte Methode wurde von Günther [1] skizziert: Da die 
Eigenwerte von (U(z)) oi,... , sind, wird die p-te Potenzsumme der Wurzeln in 
»(®) = 0 gleich der Spur!) von (U(x))” sein, oder, da sie unabhängig von x ist, gleich 
der Spur des konstanten Gliedes in der Laurentreihe für (U(x))’. Dies kann nach 
dem vorhergehenden leicht angegeben werden, weil U(x) eindeutig ist und daraus 


(30, 23) (U(z))” = Y(z, ze”) 
folgt. Hiernach können die Koeffizienten in dem Polynom o#(w) rational gebildet werden. 
Da wir nun Mittel besitzen, um die Laurentreihen für die eindeutigen Funktionen, 


die in Y,(x) und Y,;(x) eingehen, zu bilden, können wir durch den im letzten Para- 
graphen geschilderten Prozeß auch die Laurentreihe für V(x) bilden. 


$ 31. Unwesentliche Singularitäten. Die kanonische Gleichung. 


In dem speziellen Fall, daß P(x) eine isolierte Singularität in einem einfach zu- 
sammenhängenden Bereich X hat, kann man Genaueres über die Fundamentallösung 
erfahren. Liegt die Singularität in x = 0, so wird die Lösung nach dem vorhergehendem 
von der Form z4V(x) sein, wo V(x) eindeutig in Yist; weiter ist sie regulär, ausgenommen 
in dem Punkte 0, wo sie im allgemeinen eine wesentliche Singularität hat. Doch kann 
es zuweilen vorkommen, daß sie auch in 0 regulär ist, und man sagt dann, daß die Gleichung 
hier eine unwesentliche Singularität ?) hat. 


1) Unter Spur einer Matrix versteht man die Summe ihrer Diagonalelemente. Diese ist gleich der Summe der 


nach der Multiplizität aufeezählten Eigenwerte. 
2) Bieberbach, L.: Theorie der Differentialgleichungen, 2. Aufl. Berlin 1926. S. 187. 
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Trifft dies zu, so zeigt (18, 11) und (8,6), daß 
P= D,Y = Q(D.x4 | V") = a(“ | v") 
oder - 
(31,4) zP=V"AV+avv' 


ist, und hieraus erkennt man, daß, wenn V(0) invertierbar ist, die Gleichung kanonisch!) 
in x = (0 sein wird, d.h. die Koeffizienten haben in diesem Punkte einen einfachen Pol. 
Ist V(0) nicht invertierbar, so besteht eine lineare Beziehung zwischen ihren Spalten 


29,,(0)4; = 0, b am 1, ...h. 
k sei so gewählt, daß A, == O0 ist, und A bezeichne die Matrix, die dadurch entsteht, daß 


die Elemente in der k-ten Spalte der Einheitsmatrix durch #,,... 4, ersetzt werden. 
Dann wird die k-te Spalte in V(z)A offenbar einen Nullpunkt in O0 haben. Folglich ist 
(31, 2) Vv= Van, 


wenn wir mit E; = (e?) die Matrix bezeichnen, in der nur das Element e — 4 von 


0 verschieden ist. Nun ist Det. A = 4, #0, so daß 
(31,3) Det. V, = "# Det. y 


in = 0 eine Nullstelle von niedrigerer Ordnung als Det. V hat und eventuell #0 ist. 
Wiederholte Anwendung dieser Schlußfolgerung zeigt, daß man V als V,R schreiben 
kann, wobei R ein Polynom und V, invertierbar in 0 ist ?), so daß Y, = x4V, eine ka- 
nonische Gleichung befriedigt. 
Nun sagt man, daß zwei Gleichungen von derselben Art ?) sind, wenn die eine aus 
der anderen durch eine Transformation von der Form Y = Y,R hervorgeht, wobei R 
ein Polynom ist. Wir haben damit gezeigt, daß eine Gleichung mit einer unwesentlichen 
Singularität in x = O von derselben Art ist wie eine ın diesem Punkte kanonische Gleichung*). 
Andererseits: Hat P(x) einen einfachen Pol in O0, so wird für hinreichend kleines & 
|xP(x)| <cE 
sein, wobei c ein passend gewählter Skalar ist, und hieraus folgt (vgl. (25, 7)) 


Et 
| 


IYla, 2)| <1+ er _gE 


Wenn x — 0, kann also Y(a, x) nicht stärker wachsen als eine Potenz von x, so daß V(x) 
keine wesentliche Singularität in O0 haben kann. Durch passende Wahl von A ist V(r) 
also regulär in 0, d.h. die Gleichung hat hier eine unwesentliche Singularität °). 

Um eine Fundamentallösung unter der Form xz4V(x) zu finden, untersuchen wir 
zuerst die Windungsfunktion. Setzt man 


(31, 4) Pla) = + Fa), 





1) Schlesinger [1], S. 164. 
2) Vgl. Birkhoff [2], S. 388. Ein etwas umständlicherer Reduktionsprozeß, der im Prinzip von Kronecker [1], 
5. 385386 herrührt, gibt ein schärferes Resultat: V=|IV, x’ R, wo V,(0) invertierbar, R ein Polynom mit einer 
konstanten Determinante + 0 und 4 eine positive ganzzahlige Diagonalmatrix ist (vgl. Schlesinger [1], S. 150—154). 
3) Schlesinger [1], S. 105. 
*) Dieser Satz ist aufgestellt von Sauvage [2] und Koenigsberger [1], S. 449, aber zuerst vollständig be- 
wiesen von Horn [2]; einen einfacheren Beweis gab Schlesinger (vgl. Fußnote 3). 
5) Sauvage [1]. — Dieser Beweis rührt im wesentlichen von Fuhr [2] her. Früher schon gab Birkhoff [1] einen 


fast ebenso einfachen Beweis. 
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wobei F(x) regulär in 0 ist, so kann man schreiben 
y 


(31,5) U(x) = fa 1 & u. zF(at)) aı) ; 


U(x) ist offenbar auch regulär in O0, und es ist 

(31, 6) U(0) = dm, 
Nun ist U(x) = U(a), wenn weder a noch x singuläre Punkte für die Gleichung sind. 
Hieraus darf man nicht mit Volterra [2], S. 26 schließen, daß U(z) — U(0), aber jeden- 
falls ist 

(31,7) Det. (U(x) — w) = Det. (e"4 — »), 
so daß die Eigenwerte o,,...,o, von U(x) = 510g U(x) sich von den Eigenwerten 
von A nur um ganze Zahlen unterscheiden. 

Wenn #o hinreichend groß ist, kann der Integrationsweg im Ausdruck (29, 4) für 
Y(x; o) zu den Strecken (x, 0) und (0, x’) zusammengezogen werden. In das erste Inte- 
gral führen wir ein, daß 

Y(t, x) = Y(t, z)U(x) 
ist, und im anderen ersetzen wir it durch t’ = - e?*‘, wozu wir bemerken, daß 
re, 2)= Fl) 

ist. Damit bekommen wir 


zT 


(31, 8) Y(z; 0) = ft" Yt, a)dt. 
0 
Nun ersetzen wir £ durch 22 und nehmen im Produktintegral Y(zt, x) eine 2-Trans- 


formation mit t* vor; das ergibt 
1 2 
(31, 9) Y(z; 0) = ae [ #4 / (1 + 274 Fler) 7-4 dı)dt. 
’ i 


Setzen wir hierin die Peanosche Reihe für das Produktintegral ein, so finden wir folgende 
Reihenentwicklung 


Y(2;0)= D/ S.(2; 0); 
=(0) 
(31,10) | u+i) $ 
S,(x; 0) = art [ 9" Flan)u" 1%... Fler) di,...dıud. 
0 





\ 


Die (u + 1)-fache Integration erstreckt sich über den Bereich 
(31, 11) 1-3 5% 


und kann durch u + 1 iterierte Integrationen über die Intervalle 0 <ti<r,0<r, Sr, 
.:„0SsSrsi1 ersetzt werden. Führt man die erste aus und ersetzt danach r, durch 


T) 77T, durch 7, - 7, usw., so bekommt man 
1 1 
(31, 12) S.(2; 0) = "oe — A) f -.-[ Flar...7)1i  Flar... er, ,. 
0 0 
Far) Er Adır, ...dr,. 
Wenn man nun das Residuum von S,(z; o) in einem Pol bilden soll, kann man 
es z. B. in eine Potenzreihe um O0 entwickeln !). Ist 





& 
!) Es sei bemerkt, daß C,(o) = „2. @) im wesentlichen die explizite Lösung der Rekursionsformeln bei 


Horn [1], S. 395, (5) darstellt. 
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(31, 13) F(2) =2 A,a”, 
so wird 
(31, 14) Sulz; 0) = 2, (o)a’*”, 
wobei 
1 1 
= ea". SA 
at tr,=n 0 0 
(31,15) / DT a 
= I (e-A) "Aue + mA) "Aue. An let at mA)”. 
' tee tryma 


Zu einem anderen anwendbaren Ausdruck kommt man durch Deformation jedes 
der w Integrationswege zu einer Schleife von 1 um 0; setzt man 


(31, 16) F,(x) = F(x) EN 1)! 
so kann man schreiben: 
(31,17) Sa; 0) = at"(@— A) f... [Fılan...)u Alan... )2t 
1,(0),1 1,(0),1 
...Fılat,) "a 1. FR 


Hiermit haben wir S,(x2; o) für alle x in dem Regularitätsbereich von F(x) in der ganzen 
o-Ebene analytisch fortgesetzt und können nun die Partikularintegrale Y,;(x) bilden. 


$ 32. Die normale Gleichung. 

Wir nennen eine Gleichung normal in x = (0, wenn sie eine Fundamentallösung 
von der Form zP - V(x) hat, wobei V(0) eine invertierbare und P eine normale Matrix 
ist, deren Eigenwerte keine ganzzahligen Differenzen + 0 bilden. 

Man erkennt leicht, daß eine Gleichung mit einer unwesentlichen Singularität in 
0 von derselben Art wie eine in O normale Gleichung ist. Wenn die Eigenwerte von P 
ganzzahlige Differenzen + 0 bilden, kann man nämlich diese Matrix in P, + 4 teilen, 
wobei die Diagonalelemente in P, keine ganzzahligen Differenzen + 0 bilden und 4 eine 
positive ganzzahlige Diagonalmatrix ist. Aber da P, und A permutabel sein müssen, ist 


5 Er DB 
und den Faktor x° können wir vor der Konstruktion von AR(x) in V(z) einziehen. 
Ob eine gegebene Gleichung normal ist, erkennt man nach dem folgenden Satze: 


I. Eine kanonische Gleichung ist normal in 0, wenn die Eigenwerte vom Residuum 
A des Koeffizienten keine ganzzahligen Differenzen +0 bilden. Ist diese Bedingung er- 
füllt, so hat die Gleichung eine Fundamentallösung von der Form x V(x), wobei V(0) =1 ıst. 


Daß diese Bedingung notwendig ist, sieht man daraus, daß nach (31, 1) 
(32, 1) A = Res P(xz) = V'(0) PV(0) 
z=0 


ist. Um zu zeigen, daß sie auch hinreichend ist, betrachten wir das Produktintegral 


f 


V(x;&) 


l 


| (A + rt Füi)ıT" dı) x“ 
(32, 2) { . 


mM 


= ja = Pe a A) ı" dt), 
0 





| 


Journal für Mathematik. Bd. 171. Heft 2, 15 





114 Rasch, Zur Theorie und Anwendung des Produktintegrals. 


wobei & einen komplexen Parameter mit hinreichend großer reeller Komponente be- 
zeichnet (vgl. 14, I). Für die Glieder in dessen Peanoscher Reihe 


(32, 3) V(2;8)=2T,(2; 8) 
gilt infolge (19,9) die Rekursionsformel 


1 
(32, 4) Ta; &) = a". (re Tar;e)Fla)i"d, ui. 
0 


Wird 
(32, 5) Ta; 8) = #etVG,(a; $) 
und 
(32, 6) F.(2; 5) = Gu-1(2; &) F(x) 
gesetzt, so kann man schreiben: 
1 
(32, 7) RUHE IEE Kung 7977 37 ee 2 
0 
Ausgehend von der Potenzreihe 
1) k—1 
(32, 8) F(x) = Ss Aur=SAn+ x* R,(z) 


suchen wir die entsprechenden Reihen 
Rn ri 5 Due 
(32, 9) Fta;!)= 2A (= SA (dr reR (2) 
und 
BEE 2" SB ar ee 
(22,10) Wk) edrrrd lat). 


Die Koeffizienten werden rekursiv durch die Formeln 


(32, 11) ANE) = SB E)Ar-ı 
und 
1 
(32, 12) Be) = JE 


bestimmt. Unter Anwendung der dritten kanonischen Form für {4 kann die zweite 


dieser Formeln 
i —y—1 An A” (E) A, 
(32,13) BY) = : 
z\ 1 a ee er 


geschrieben werden, wobei die Summation sich über alle vorkommenden Werte von 
h,j, vund Aerstreckt. Die Koeffizienten sind also rationale Funktionen von £ mit Polen in 





0, 0,7% 


(32,14) = I  —l himhm rl. Behr... 


Weiterhin ist das Integral 


1 
0 


konvergent für 























EEE te 
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(32, 16) RE> En o = Max Re, — o,), 
u (h, ;) 
und hieraus schließt man durch Induktion, daß die Funktionen F,(x;&) und G,(x; €) 
meromorphe Funktionen von & mit den Polen (32, 14) sind, falls x im Inneren des prinzi- 


palen Mittag-Lefilerschen Sternes © von F(x) liegt. 


Nach der Annahme, die wir über die Eigenwerte von A gemacht haben, kann der 
Punkt & = 0 nicht unter diesen Polen vorkommen, und er kann auch kein Häufungs- 
punkt für diese sein. Folglich sind sämtliche Glieder in der Peanoschen Reihe (32, 3) 
in einer gewissen Umgebung von & = 0 regulär. Weiterhin werden wir zeigen, daß die 
Reihe dort gleichmäßig konvergent ist, so daß auch V(x; &) in & = 0 regulär ist, wenn 
x im Innern von © liegt und arg x beschränkt ist. 


Aus der Halbebene RE > — 1 + n > — 1 werden eventuelle Pole (32, 14) heraus- 
geschnitten, und die ganze Zahl qg wird > 7 @ + 2e), gewählt, wo e eine positive Zahl 


ist, die wir gleich besprechen werden. Wir wollen nun zwei Folgen von positiven skalaren 
Funktionen (f,(2;&)) und (g,(x;&)) konstruieren, die für « =g die Ungleichungen 


(32, 17) FAST HE, IC, Sg, ee 


befriedigen, wenn z in einem geschlossenen Bereich in © und £ in einem geschlossenen 
Bereich in der präparierten Halbebene liegt. Aus E?= E und (32, 6) folgt, daß wir 


(32, 18) 9) = g,_.l; fe) 
setzen können, indem 
(32, 19) |F(z)| Ss f(@)E. 


Aus der dritten kanonischen Form für t4 erkennt man, daß, wenn po’ und o’’ die Eigenwerte 
von A mit dem kleinsten bzw. größten reellen Teil sind, man zu einer beliebigen positiven 
Zahl e einen konstanten Skalar x wählen kann, für den 


ısVa®E, tl syar® E 
ist, mit 0 <z?<4. Mit dem gewählten g ist weiterhin 


Daraus folgt, daß wir für vu >21 
1 
(32, 20) BRUR C51 3 El Ran FFIR C 723. 
0 


wählen können, indem wir der Kürze halber R& = £, setzen. Durch Induktion ist leicht 
zu zeigen, daß dies 


(& rafler)dr)e”" 


(32, 21) +85 8) > | u =: I f(at)g (at; E) dt 





ergibt. Da F(x) und G,(x; £) in den betrachteten Bereichen in der x- und &-Ebene regulär 
sind, können f(x) und g,(2; &) gleichmäßig beschränkt gewählt werden, und damit er- 
halten wir eine gleichmäßig konvergierende Majorantenreihe für die Peanosche Reihe 
(32, 3). 

Für 

(32, 22) Y(z; 8) = 2*V(2; 8) 


15* 
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gilt infolge (32, 2) und (18,11) die Relation 
D.Y(z; &) = Q(DAx4V (x) 24) | 24) = QxttFlx) 2-4] 2-4) 


BE - z + #F(e). 
Da V(x;&) für die angegebenen x in &= 0 regulär ist, wird 
(32, 24) Y(z) = Y(2;0) = =“ V(2) 
die Gleichung 
(32, 25) Y'(xz) = Y(x) P(x) 


befriedigen. In der Umgebung des Punktes x = 0 wird die Reihe 
(32, 26) V(a) = 2 T,(2;0) 

gleichmäßig konvergieren, so daß V(x) in O0 regulär ist und 
(32, 27) v0), =1. 

Damit ist der Satz I bewiesen !). 


Wir sind hiernach imstande, eine gegebene kanonische Gleichung, deren Koeffizient 
das Residuum A hat, in eine normale Gleichung zu transformieren. Zuerst wird eine in- 
vertierbare Matrix A gesucht, die A in seine Normalform P transformiert; wird Y durch 


YA ersetzt, nimmt die Gleichung die Form 
(32, 28) I!’ = H- 4 F) = YP, 
x 


an, wo F(x) regulär in O ist. Weicht nun einer der Eigenwerte von A, z.B. 0, von einem 
anderen um eine ganze positive Zahl ab, so teilen wir P in eine Relativmatrix 


(32, 29) P = Ps 5) 
0 P, 
ein, woP, =, + ©, und die Eigenwerte von P, 0a, 0, sämtlich # o, sind. Wird 
die Gleichung mit der Relativmatrix 
x 0 
(32, 30) x=(b ı) 


transformiert, wobei x ein Skalar von gleicher Ordnung wie P, ist, bekommt man eine 
Gleichung mit dem Koeffizienten 


| 
(32, 31) P,= PX) = a X) + XFX”" 
Wird F auf gleiche Weise wie P in 
Fu F 
32, 32 pi zu m) 
F3ı Foo 


!) Im übrigen sei bemerkt, daß der Beweis sich leicht so modifizieren läßt, daß die Regularität von F(x) in 0 
dadurch ersetzt wird, daß F(x) für x — 0 längs eines gegebenen Radiusvektors oder in einem gegebenen Winkelraum 
asymptotisch durch eine Potenzreihe 


[+] 
Fix) 3 Au“ t?=!, 5>0 
x=(0 


dargestellt werden kann. Man findet dann V(x) —1 für x — 0 und allgemein 


V(z) Le 0) -+ O(2”°). 
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eingeteilt, so wird 


ö Fı #Fj 
32, 33 XFX =|1 
\ —Fa Fa 
und da X permutabel mit P ist, erhält man 
Fr )= 7 - 
(32, 34) a(- %)=-—(', 5)" 
Die neue Gleichung ist also auch kanonisch und das Residuum von P, 
- P, u 1 0 
EB = (Fo) PR} 


hat die Eigenwerte o, —1,o,...0o,. Wenn sie keine ganzzahligen Differenzen haben, 
dann ist die Gleichung normal, und sonst behandelt man sie genau wie die ursprüngliche. 
Man sieht, daß man mit einer endlichen Anzahl von Transformationen zu einer normalen 


Gleichung kommen kann. 
Wir haben oben ($ 31) bemerkt, daß U(x) zuweilen nicht mit U(0) äquivalent ist. 
Dieser Fall kann nicht eintreten, wenn die Gleichung normal ist. Dann ist nämlich 


2niA 


(32, 36) U(x) = V'(2) &""V(a) me“, 
wo A das Residuum des Koeffizienten und V(0) = 1 ist. Wenn sie nicht normal ist, hat 
sie eine Fundamentallösung von der Form 

(32, 37) Y(z) = 24 V(x) R(x), 
wobei R(x) ein Polynom ist, V(0) = 1 und A das Residuum des Koeffizienten in einer 
normalen Gleichung von derselben Art wie der gegebenen bezeichnet. Folglich ist 


(32, 38) U(z) = R(&)V (x) E** V(x) Rx) mei, 
während 
(32, 39) U(0) = edri4, 


und obgleich e*"* und e’”“ dieselben Eigenwerte haben, brauchen sie keineswegs äqui- 
valent zu sein. Wenn man den regulären Teil #(x) des Koeffizienten nicht mit in die Dis- 
kussion beziehen will, läßt sich nur folgendes behaupten: Die Normalform von A, 


Pr +.:0 
(32, 40) -( Zraees 


sei so geschrieben, daß 1. die P, je nur einen Eigenwert o, haben, 2. 0, + oe, für j + 


ist, und 3. in der Anordnung P,,P,,..., Pm diejenigen P; unmittelbar aufeinander 
folgen, deren Eigenwerte ganzzahlige Differenzen besitzen, und zwar nach der Größe 


der Realteile der Eigenwerte geordnet. Dann wird die Normalform von A von der 
Form 


P 0 0 
(32, 41) Palit 0 9 


u. u 


sein, wobei P;,— P; ein ganzer Skalar > 0 ist und alle Elemente in T; gleich 0 sind, aus- 
genommen dem in der oberen rechten Ecke, das gleich 1 sein kann, wenn 0,=0,;,,1 ist. 
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In der Theorie über skalare Gleichungen von n-ter Ordnung behandeln Fuchs [1] 
und Frobenius [1] die Frage, wieweit ‚erwartete‘ Logarithmenpotenzen in den Lösungen 
wirklich auftreten und speziell, ob Logarithmen ganz fortfallen. Von unserer Seite ge- 
sehen liegt der Grund zum Entstehen dieser Frage darin, daß die Gleichung auf Matrix- 
form gebracht wohl kanonisch, aber im allgemeinen nicht normal ist, indem die Wurzeln 
einer gewissen algebraischen Gleichung, die unserer charakteristischen Gleichung 


Det. (A— o) =0 
entspricht, ganzzahlige Differenzen haben können. Die aufgestellte Bedingung für die 


Verteilung der Logarithmenpotenzen gibt in Wirklichkeit eine Transformation an, die 
die kanonische Gleichung in eine normale überführt. 
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Ein „belastetes“ Variationsproblem. 


Von H. Boerner in München. | 


Wenn man kompliziertere Variationsprobleme behandelt, z. B. das Lagrangesche 
Problem, dann fügt man bisweilen zum Integral noch eine Funktion der — auf irgend- 
welchen Mannigfaltigkeiten beweglichen — Endpunkte hinzu !); man erhält dann ein 
Zusatzglied in den Randbedingungen, dessen Mitnahme die Durchführung der Aufgabe 
nicht weiter erschwert. Interessanter werden die Verhältnisse, wenn man eine solche 
„Last“ irgendwo im Innern anbringt, was soweit ich sehe bisher noch nicht geschehen 
ist. Im folgenden soll kurz das allereinfachste Variationsproblem in dieser Weise behan- 
delt werden und zwar, weil es sich dabei um eine „Ecke‘‘ handeln wird, in Parameter- 
darstellung. Die Anregung zu einer Problemstellung dieser Art verdanke ich Herrn Pro- 
fessor Lichtenstein }. Die Übertragung auf Doppelintegrale dürfte nicht schwer 
sein; es handelt sich dann um das Extrem der Summe eines Doppel- und eines einfachen | 
Integrals, ein Problem, das Hilbert einmal behandelt hat ?), aber wiederum nur für den } 
Fall, daß das einfache Integral über einen Teil des Randes erstreckt wird. ; 





1. Es sei T ein einfach zusammenhängendes Gebiet der x, y-Ebene, L ein stetig | 
gekrümmter Querschnitt, der 7 in zwei Teilgebiete 7, und 7, zerlegt. Ist t die von 
irgendeinem Anfangspunkt aus gezählte Bogenlänge auf ZL, so seien x = £(t), y = nt) 
die Gleichungen von ZL; £ und n sind mit ihrer ersten und zweiten Ableitung stetig; wir 
setzen noch &’(t) = cost, n’(t) = sin r; die positive Richtung soll so gewählt sein, daß 
T, zur Linken, T', zur Rechten liegt. Ferner sei auf Z eine zweimal stetig differenzierbare 
Funktion ®(t) gegeben. 


Die Funktion F(z, y; x’, y’) ist in 7 und für alle x’, y’ definiert und genügt den 
üblichen Differenzierbarkeits- und Homogeneitätsbedingungen °). Es ist 








(1.1) x Fr +yFy=F 
und 
Fr Fr ’ F 4 ‚ : 
(1.2) Du = F,(2, y; ®, y'); 


die E-Funktion ist durch 
(13) Eia,y;20',y';2,y)= {Fr(z, y; ®, y) — Fx(z, y; ©, y’)} © 
+ {Fyla, y; 2’, y') — Fa, y; &, y’)} y' 
definiert. 


Das Variationsproblem lautet so: Unter allen Kurven C:z=xf(s), y= y(s) 
(s irgendein Parameter) mit abteilungsweise stetiger Tangente, die einen festen Punkt P, 





ı) 2.B.O. Bolza, Math. Ann. 74 (1913), p. 430; M. Morse, Am. Journ. of Math. 58 (1931), p. 517. 
2) Math. Ann. 62 (1906), p. 365. 
®) Vgl.z. B.O. Bolza, Vorlesungen über Variationsrechnung, Leipzig 1909, p. 193#f. 
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in 7, mit einem festen Punkt P, in 7, verbinden und den Querschnitt Z nur einmal 
treffen und nicht berühren, diejenige zu finden, die der Summe 


(1.4) Je = [Fz, y, #', y') ds + ®(t) 


den kleinsten Wert erteilt. Hierbei ist i der Parameterwert auf L, der zum Treffpunkt 
gehört. 

Wie üblich, betrachten wir eine Kurve C und stellen notwendige und hinreichende 
Bedingungen dafür auf, daß sie die Lösung liefert. Der Parameter s soll die Bogenlänge 
bedeuten und läuft von s, bis s,; wir setzen x’(s) = cos ®, y’(s) = sin 9. Der Treffpunkt 
mit L sei P, mit den Koordinaten x,, 4, und den Parameterwerten s, auf C und t, auf Z. 
Da C ım allgemeinen im Punkt ?P, eine Ecke haben wird, führen wir für die beiden Teil- 
bögen von P, bis P, und von P, bis P, die Bezeichnungen C',, bzw. C'y, ein und nennen 


die Winkel, die ihre Tangenten in P, mit der x-Achse einschließen, 9, bzw. 9,; 7, ist der 
Winkel zwischen der x-Achse und der Tangente von ZL in P,. 


Betrachtet man zunächst Vergleichskurven, die Z ebenfalls in ?, trefien, so erhält 
man zwei gewöhnliche Variationsprobleme mit den festen Endpunkten P,, P, bzw. 
Po Pa; die Bögen C,, und C,, müssen also der zu F(x, y; x’, y’). gehörigen Eulerschen 
Differentialgleichung genügen, und wir nehmen an, daß die bekannten Bedingungen 
von Legendre, Jacobi und Weierstraß in der folgenden Weise erfüllt sind: 

(1.5) es ist FR, > 0 im ganzen Intervall (s,, s2); 

(1.6) dieses abgeschlossene Intervall enthält keinen zu ?, konjugierten Punkt; 

(1.7) die E-Funktion, gebildet für die Richtung # und irgendeine andere Richtung, 

soll nie negativ sein und höchstens im Punkt P, außerordentlich verschwinden; 

(1.7) wird nur verlangt, wenn es sich um starkes Minimum handelt. Ferner soll 
die Extremale C außer in P, keine Ecke besitzen, eine Voraussetzung, die ohnehin üblich 
ist, wenn man von schwachem Minimum redet; im Falle des starken folgt sie schon 
aus (1.7). 

Für die Funktionswerte im Punkte ?, führen wir noch die abkürzende Bezeichnung 
F(x,, %0; c08 9, sin 9) = F(9) ein; entsprechend für die partiellen Ableitungen von F 
und für F,; die E-Funktion heiße 

(1.8) E(9, 9) = E(x,, Yo; 608 ®, sin 9; cos d, sin 9). 

Betrachtet man jetzt beliebige Vergleichskurven, so handelt es sich um die Summe 
von zwei Integralen mit gemeinsamem auf Z beweglichen Endpunkt, vermehrt um ®. 
Wir werden darum den Parameter so wählen, daß er ebenfalls von s, bis s, läuft und daß 
der Wert s, zum Treffpunkt gehört. Jede Vergleichskurve, deren Treffpunkt von ?, 
verschieden ist, kann als Mitglied einer Kurvenschar 

(1.9) = pls), y=ylsl) 
dargestellt werden, so daß 
ps, to) = 2(s),  Yls to) = Yls) 

Ist) = lt), Yon) = ne). 
Für s = s, ist also (die Indizes bezeichnen die partiellen Ableitungen) 


(1.10) 


%,=c08T, y=Ssinrt. 


Die Eckenbedingung erhält man, indem man die erste Variation bildet, d.h. die 
Summe J, als Funktion von t betrachtet, differenziert, i = t, setzt und die Tatsache 
benützt, daß die Kurve C aus Extremalenbögen besteht, und sodann die erste Variation 


gleich Null setzt. Man erhält ohne weiteres 
Journal für Mathematik. Bd. 171. Heft 2. Ib 
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(1.11) {Fztdg) — Fr(do)} 608 To + {Fy(d) — Fy9o)} sin Tg = Dt) 

2. Um die Gleichung (1.11) zu diskutieren, bezeichnen wir mit A irgendeine Zahl 
und betrachten die Gleichung 

(2.1) {F(&) — Fr’(a)} cos ty + {F,(&) — Fy(a)} sin zu = A 
für irgendwelche Winkel «, x, die zwischen 7, und 7„— x liegen. Für die partielle Ab- 
leitung der linken Seite von (2.1) nach & erhält man mit Berücksichtigung von (1.2) 


den Ausdruck 
— F,(a) sin (n— a), 


entsprechend 

+ F,(&) sin (— &) 
für die Ableitung nach a. Erteilt man also dem Winkel «& einen festen Wert, für den 
F,(&) > 0 ist, und beschränkt den Winkel & auf diejenigen Intervalle, in denen F,(&) > 0 
gilt, so wird in jedem dieser Intervalle die linke Seite von (2.1) monoton mit & wachsen. 
Ist A eine geeignete dem Wertebereich der linken Seite von (2.1) angehörige Zahl, so gibt 
es eine Anzahl Lösungen & von (2.1), darunter, wenn nur |A | genügend klein ist, gewiß 
eine ausgezeichnete, bei der x in demselben Intervall liegt wie x. Variiert man A, so 
ändert sich & im gleichen Sinn; wird dabei A = (0), so ergibt die ausgezeichnete Lösung 
% = x, also eine kontinuierliche Extremale des gewöhnlichen Variationsproblems. Aus 
den andern Lösungen wird dagegen im allgemeinen keine geknickte Extremale des ge- 
wöhnlichen Problems; denn die Erdmannschen Bedingungen 

(2.2) Fy(&) = Fr(a), Fyla) = Fyla) 

verlangen ja viel mehr als (2.1) mit A = 0, und daher kommt es, daß die Lösungen 
dieser Gleichungen stets isoliert liegen, während man aus jeder Lösung von (2.1) durch 
variieren von & beliebig viele neue erhalten kann. Es gilt aber der Satz: 


Ist für die Winkel x, & die Gleichung (2.1) mit A=0 erfüllt und sind die durch & 
und x ın T, bzw. T, bestimmten Extremalenstücke in der Umgebung von P, stark, so gilt (2.2). 


Zum Beweis betrachten wir den Ausdruck 
y(ß) = {Fr(&) — Fr(a)} cos P + {Fr(&) — Fr(a)} sin P 


für alle Richtungen f. y(ß) verschwindet nach Voraussetzung für? = „und? = n—ın 
und ist, falls (2.2) nicht erfüllt ist, im Intervall , — x < 8 < r, entweder immer positiv 
oder immer negativ, die Zahlen y(x&) und y(%&) sind also beide positiv oder beide negativ. 
Nach (1.3) ist aber mit der Schreibweise (1.8) 


y(a)=E(6,0), yla)=—Ela,e). 


Wir haben jedoch 
E(%,&)>0, E(&,a)>0 


vorausgesetzt, also muß y(x) = y(a) = 0 sein, woraus (2.2) folgt. 

Ist insbesondere P, ein regulärer Punkt !) des gewöhnlichen Variationsproblems, 
also F, > 0 für alle Richtungen, dann folgt aus (2.1) mit A = 0 stets nicht nur (2.2), 
sondern sogar & = a. 

3. Um die Weierstraßsche hinreichende Bedingung abzuleiten, braucht man ein 
Feld von geknickten Extremalen. Wegen (1.6) kann man ein Feld von kontinuierlichen 
Extremalen konstruieren, das in 7, durch die Gleichungen (1.9) gegeben ist, den Be- 


ı) Nach der Definition von Hilbert, s. z. B. Bolza, Vorlesungen p. 214. 
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ziehungen (1.10) genügt und den Punkt P, sowie ein Stück der Kurve Z in der Umgebung 
von P, im Innern enthält. Das Feld ordnet jedem Punkt des von ihm bedeckten Ge- 
bietes ein Gefälle zu, das wir mit ®(s, t) bezeichnen können. Die Gleichung 


(3.1) {Fy(z, y; cos d, sin 9) — Fy-(x, y; cos ®, sin O)} cos r 
+ {F, (x, y; cos d, sin 9) — F,(x, y; cos ®, sın Ö)} sin r = ®’(t)!), 
in der x = £(t), y = nlt), d = %(s, t) zu setzen ist, wird für i = t, durch den Winkel 9, 
befriedigt, und in diesem Punkt ist, wie wir am Anfang des vorigen Abschnitts sahen, 
die Ableitung nach # von Null verschieden. Also kann man nach # auflösen und in den 


so erhaltenen Richtungen, die wir mit #(s,, t) bezeichnen, von den Punkten von L aus 
ins Gebiet 7, hinein Extremalenstücke konstruieren, die unser geknicktes Extremalen- 
feld (1.9), (1.10) vervollständigen. Wir beschränken dabei t auf eine so kleine Umgebung 
von t,, daß im ganzen Feld die Legendresche Bedingung (1.5) erfüllt bleibt. Nehmen 
wir noch an, daß die Jacobische Bedingung für das „belastete‘‘ Problem, die wir im 
nächsten Abschnitt ableiten werden, erfüllt ist, so ist es möglich, das Feld so anzulegen, 
daß es auch die Kurve C,,, also die ganze Extremale C im Innern enthält. 


Ist in diesem Felde überall 
E(x, y; cos ®, sin 9; cos ®, sin d)>0, 
wenn man 9 = Ds, t) und für d irgendeine Richtung + 9 wählt, so liefert die Extremale C 
ein starkes Minimum. 

Zum Beweis betrachten wir irgendeine Kurve C:x=ä(o), y= (ce) mit ab- 
teilungsweise stetiger Tangente, die ?, mit P, verbindet, ganz im Felde verläuft und die 
Kurve Z im Punkt P, mit dem Parameterwert t trifft. Das „Hilbert-Integral“ 
[tFr(&, y, cos d, sin 9) £’ + F,’(z, y; cos ®, sin 9) y’} de 
f 


(d ist das Gefälle des Feldes im Punkt &, 7), das sich für C == C auf die Größe Je — ©(t,) 
reduziert, ist nur in jedem der beiden Teilfelder links und rechts von Z vom Wege unab- 
hängig. Aber der Integrationsweg darf jeweils auch Teile von Z enthalten. Integriert 


man nun über Z vom t bis t, mit den Grenzwerten des Gefälles von der Seite 7, und dann 
zurück von t, nach t mit den überstrichenen Werten, so erhält man wegen (3.1) einfach 


(ou dt = Bft) — ©(t,). 


Daraus folgt, daß die Summe 
J{Fr(&, 7;cos d, sin BP) E' - F(&, 7; cos d, sin PB) y’}do | ®ft), 


die für C = C in J« übergeht, überall vom Weg unabhängig ist. Diese Summe ist also 
gleich J_, wenn C die Punkte P, und P, verbindet. Subtrahiert man sie von J;, so fällt ® 


ganz weg und man erhält unmittelbar wegen (1.1) und (1.3) und unter Benutzung 
der Homogeneitätseigenschaft der ersten Ableitungen von F die Weierstraßsche Formel 
für unser Problem: 


J.— J,= [E(&, 7; 008 d; sin 9; 8,7’) do, 
rR 


aus der wegen (1.7) die zu beweisende Behauptung unmittelbar folgt. 


ı, Wir können cos # schreiben statt y,(a,, £) usw., weil die Ableitungen F,’ und Fy’ in diesen Argumenten 
homogen vom Grad 0 sind. 
16* 
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4. Wir betrachten die Gesamtheit aller Extremalenfelder von der vorhin betrach- 
teten Art, die man in der Umgebung von P, konstruieren kann. Wenn man annimmt, 
daß auf den Bögen C,, und C,, (oder auf ihrer Verlängerung über ?, bzw. P, hinaus) 
die zu P, konjugierten Punkte P5 und P% existieren, so besitzt jedes dieser Felder auf 
jedem dieser Bögen P&P, und P,P, genau einen Brennpunkt, Q, und Q,. Wir werden 
zeigen, daß die hierdurch definierte Zuordnung zwischen den Punkten von P)P, und 
denen von P,P% umkehrbar eindeutig ist, d.h. daß wenn zwei Felder denselben Brenn- 
punkt Q, besitzen, auch die entsprechenden Punkte Q, zusammenfallen und umgekehrt, 
und ferner daß sie monoton ist, d.h. daß wenn Q, von P; bis P, läuft, sich der ent- 
sprechende Punkt Q, von P, nach P% bewegt. Das berechtigt uns dann, Q, den zu Q, 
konjugierten Punkt zu nennen. 

Betrachten wir ein bestimmtes Feld, gegeben durch die Gleichungen (1.9), (1.10), 
(3.1), und beschränken wir uns zunächst auf den in T, gelegenen Teil, s < s,. Die Funktion 


(4.1) Als) = | MP —_ y,(s, t,) cos d— g,(s, t,) sin ® 
| 9% Yı = 
ist eine Lösung der zum Bogen C,, gehörigen Jacobischen Differentialgleichung des ge- 


wöhnlichen Problems, die für den zum Punkt Q, gehörigen Parameterwert verschwindet 
und durch diese Nullstelle und durch den Wert 


(4.2) A(s,) = sin T, COS 9, — C08 T, sin 9, = sin (y— 9) , 


der von der Wahl des Feldes überhaupt nicht abhängt, eindeutig festgelegt ist. Die zu 
verschiedenen Nullstellen gehörigen Lösungen lassen sich auch durch die verschiedenen 
Werte 4’(s,) ihrer Ableitung im Punkt s, charakterisieren, die, wie man aus den Sturm- 
schen Separationssätzen leicht folgert, monoton von — © bis + © wächst, wenn 
der Punkt Q, von P% bis P, variiert. 


Entsprechendes gilt für die andere Feldhälfte. Hier ist 

(4.2’) A(s,) = sin (un — 9,), 
und der zugehörige Wert A’(s,) wächst ebenfalls von — & bis + ®, wenn der Punkt Q, 
den Bogen P,P% durchläuft. Wir werden nun eine Beziehung der Form 

(4.3) As) = AA ls) +B, A>O, 
herleiten, die allein aus dem Wert von A’(s,) den zugehörigen Wert von 4’(s,) zu be- 
rechnen gestattet, und damit werden wir bewiesen haben, was wir beweisen wollten. 


Zu dem Zweck betrachten wir die zu einem Feld gehörige Funktion #(s, t), zunächst 
wieder in 7T,. Es ist 








Y, Y, 
(4.4) O8 Gm ale, In nn 
Vp?+ y Vp+ y 
os = 005 d v —- = Ze —— sın d _ _ Pt 
! Vp2+ y VPE + y 


da die vom Differenzieren der Wurzel herrührenden Terme sich zerstören. Für s = s,, 
t = t, hat man einfach 


0» . 
(4.5) (37), = 008 Poly) — sin Dal ya) 


Differenziert man die Relationen (4.4) nach s, multipliziert die erste mit y,, die zweite 
mit — 9,, addiert und setzt s = s,,t = t,, so erhält man 


od 
(4.6) — c08 (Ty — %)) 5), = (Py)o In Tg — (Yes)o COS Ty. 
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Nun differenzieren wir (4.1) nach s und setzen s = s,; dann erscheint rechts die 
Summe der rechten Seiten von (4.5) und (4.6). Also ist 


(4.7) 2’(s,) = &).- cos (Ty— 9,) 2). 


Zu (4.5), (4.6), (4.7) analoge Formeln, die zur andern Feldhälfte gehören, erhält man, 
indem man über alle 9 und A Querstriche setzt. 
Endlich wird (3.1) nach t differenziert und i = t, gesetzt; unter Benutzung von 
(1.1) erhält man 
dcos tr | i dsinr 
ed) — Fett} () + —Rr 0) 
B Fo, dt /gmı, 
+ {Fr(d) — Frrz(Bo)} 608? Tg + {Frry(90) + Fyraldo) — Frryl®o) — Fyrz(Bo)} cos T, sin Ty 
h En =, [09 Ki) 
+ {F (do) — Fyy(9o)} sin? To + Fı(9,) sin (Tg — 9%) ,- F,(9,) sin (7, — 9,) En ) 
= ®P'(t,). 





en von der Wahl 
ot /o 


des Feldes ab; wir können A’(s,) und A’(s,) in die letzte Gleichung einführen und alles 


In dieser Geichung und in (4.7) hängen nur die Größen I und | 
0 


außer dem Term mit A’(s,) auf die rechte Seite bringen. Dann sind die Koeffizienten 
von A’(s,) und A’(s,) beiderseits positiv, man kann durch den ersteren dividieren und 
erhält die gewünschte Beziehung (4.3). 


Jetzt können wir wieder die vorgegebenen Endpunkte ?/, und P, betrachten. 
Dann lautet die Jacobische notwendige Bedingung für ein Minimum: Der auf der geknickten 
Extremalen C zu P, konjugierte Punkt darf nicht vor dem Endpunkt P, liegen. Und diesen 
konjugierten Punkt kann man ohne Feldkonstruktion so erhalten: Man bestimmt das- 
jenige Integral der zum Bogen C,, gehörigen Jacobischen Differentialgleichung, das für 
s = s, verschwindet, verfügt über den darin enthaltenen willkürlichen Faktor so, daß 
A(s,) = sin (Ty— 9,) wird, bestimmt A(s) als Integral der zu C,, gehörigen Jacobischen 
Gleichung aus A(s,) = sin (T,— ®,) und (4.3) und sucht die nächste Nullstelle von A(s). 

Man bemerke, daß die Verhältnisse in der Umgebung unserer Ecke P, viel ein- 
facher liegen als in der Umgebung einer ‚freien‘ Ecke des gewöhnlichen Problems nach 
der von Caratheodory, Dresden und Bolza entwickelten Theorie der diskontinuierlichen 
Lösungen !). Sowohl was die Theorie der ersten Variation wie vor allem die eben ent- 
wickelte der konjugierten Punkte betrifft, besitzt der Punkt ?, viel mehr Ähnlichkeit 
mit einem gewöhnlichen Punkt als mit einer Ecke des gewöhnlichen Variationsproblems. 
Es liegt dies an folgendem. Wir haben eine Zuordnung der Integrale der Jacobischen 
Gleichung für die beiden Teilbögen betrachtet. Diese Zuordnung ist so beschaffen, 
daß niemals das eine Integral in ?, verschwinden kann, wenn das andere von Null 
verschieden ist. Gerade dieses Vorkommnis ist aber charakteristisch für die Kompli- 
zierung der Verhältnisse, die man in der Theorie der konjugierten Punkte bei den 
diskontinuierlichen Lösungen des gewöhnlichen Problems hat. 


5. Man kann das zuletzt behandelte Jacobische Kriterium auch durch ein „Eigen- 
wertkriterium‘‘ ersetzen, wenn man mit Lichtenstein ?) von der Theorie der „zweiten 
Variation‘ ausgeht. 





1) Vgl. Bolza, Vorlesungen p. 372ff., wo auch die Originalabhandlungen zitiert sind. 
2) Vgl. L. Lichtenstein, Gött. Nachr. 1919, p. 161. 
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Dazu betrachten wir eine beliebige Schar von Vergleichskurven, deren Darstellung 
wir, ein wenig allgemeiner, so ansetzen: 
a = o(s,&), . y(s, &). 
Für x = a, soll die geknickte Extremale C herauskommen und wie früher fürs =s, 
die Kurve Z; d.h. man hat eine Funktion t(&), so daß 
950, %) = Elt(a)),  Yls0, a) = nltle)) 
ist. Ferner führen wir die Funktion 


Pr Pa 
w(s) = | | 
IT| yalana 
ein; für s = s, wird wie früher der zum Bogen C,, gehörige Wert durch Überstreichen 


gekennzeichnet. Es ist also, wenn man noch t’(&,) = w, setzt, 


= Yals, %,) cos d— 9,(S, &,) sin 9 


(5.1) w(s,) = wy-sin (y— 9), ls) = w- sin (un — dB). 


Die zweite Variation ist die zweite Ableitung der Größe J, die ein Minimum sein 
soll, nach «, gebildet an der Stelle x,. Die Weierstraßsche Umformung !) liefert für sie 
den Ausdruck 


(5.2) 2J, = [ {pw — kw} ds — hu; 


dabei ist 
ps) = Fılay, #y')> 0, Kı)=— Fılay, y)); 
h = {L(9,) — L(9,)} cos? ru + 2{M(9,) — M(9,)} cos r, sin Ty 
+ {M 9) — N(8,)} sin? 79 + {Fz(9,) — Fr(9)} &’(to) 
F {F (9) —Fy(Bo)} 7") —P’lto); (L(9,) — L(x, y,cos®,,sin#,) usw.). 
Die Funktionen p und k sind wie w ım Punkt s, unstetig; p(s,) ist der Grenzwert vom 
Bogen Cya her. 

Die zweite Variation soll nicht-negativ sein für alle außer in s = s, stetigen Funk- 
tionen w(s) mit abteilungsweise stetiger Ableitung, die für s = s, und s = s, verschwinden 
und für s = s, mit passendem w, (5.1) erfüllen. Betrachten wir also (nach Bliss) das 
Problem, sie zum Minimum zu machen. Dazu muß man eine Normierung einführen, 
etwa 


(5.3) [kweds+hw=1 


(indem wir uns für den Moment auf diejenigen w beschränken, für die dieser Ausdruck 
positiv ist; es sind ohnehin die einzigen, für die J, Null oder negativ werden kann). 
Es handelt sich also um das isoperimetrische Problem, die zweite Variation oder, was 
auf dasselbe hinauskommt, das Integral 


(5.4) [pw ds 


zum Minimum zu machen unter der Nebenbedingung (5.3). Dies führt uns von selbst 
zum Studium des Eigenwertproblems der Differentialgleichung 


d ’ Ihr. 


!) S. z. B. Bolza, Vorlesungen p. 224ff., wo man auch die Erklärung der Funktionen L(zy, x’y’), M, N und 
F,(zy, x’y’) nachsehe. 
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mit der Randbedingung w(s,) = w(s;) = 0 und den Unstetigkeitsbedingungen (5.1) und 
(5.5) Pf(so) sin (Tu — 95) » w’(sg) — P(So) sin (Tg — Bu) - w’(s,) + Ahw, = 

im Punkt s=s,!). Dieses Problem kann man mit der Lichtensteinschen Methode ?) 

behandeln; diese liefert die Existenz von unendlich vielen Eigenwerten, im allgemeinen 


positiven und negativen, A, und von zugehörigen normierten Eigenfunktionen g(s), 
die den Orthogonalitätsbeziehungen 


fh op ds -- hyi pe) — Öj 2 3) 
8, N 

genügen. Zu jeder zulässigen Funktion w(s) — die Normierung (5.3) und die im Zu- 

sammenhang damit gemachte Einschränkung ?) fällt wieder fort — kann man die ‚‚Fourier- 


koeffizienten‘“‘ 


83 
x; — [kwgWds + hw, pV) 


8] 


einführen und erhält die beiden Beziehungen 


a ’ an 
[kw ds + hw= > un 
$%ı j 14;| ‚ 


und 


Spw?ds>22|2| x, 
s J 


ı 


aus denen sofort die Lichtensteinsche Formel für die zweite Variation folgt: 


1,25 I1s1(1-,.)e 


Man braucht also bloß den kleinsten positiven Eigenwert /, zu betrachten (man 
setze A, = + @, falls nur negative Eigenwerte existieren). Es muß jedenfalls }, > 1 
sein, damit die zweite Variation nicht negativ werden kann; das ist die der Jacobischen 
entsprechende Bedingung. Man zeigt leicht, daß mit der zweiten Variation auch die 
Gesamtvariation negativ wird, sobald man |x — a,| genügend klein wählt. Darüber 
hinaus ergibt sich, daß die schärfere Bedingung /, > 1 für ein schwaches Minimum 
hinreichend ist. Ist A, = 1, so müssen die höheren Variationen herangezogen werden. 





1) Für A= 1 handelt es sich um die Jacobische Differentialgleichung für die Bögen C'yo, Ca, und (5.5) ist 
dann natürlich im wesentlichen identisch mit (4.3). 

2) Vgl. L. Lichtenstein, Rend. eirc. mat. Palermo 88 (1914), p. 113. 

®) ö,, ist das Kroneckersche Symbol. 


4) Ist g(s) eine Eigenfunktion, also eine Extremale des isoperimetrischen Problems, und A der zugehörige 
Eigenwert, so ist stets 


83 2 
[p»* af kt ds + ho = 0. 
8 8 


Daraus folgt, daß A, das Minimum von (5.4) unter der Nebenbedingung (5.3) ist. Bezeichnet A_, den kleinsten 
negativen Eigenwert, so ist — A_, das Minimum von (5.4), wenn man die rechte Seite von (5.3) durch — 1 er- 
setzt. Die entsprechenden Minima der zweiten Variation sind dann A, —1bzw.|i_,' +1. 


———— nn m 


Eingegangen 6. April 1933. 
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Zu einem Satze von Hazzidakis über gewisse Systeme 
linearer Differentialgleichungen. 


Von Friedrich Dingeldey in Darmstadt. 





Im 90ten Bande (1881) dieses Journals (S. 80—82, [1880]) hat J. N. Hazzidakis 
einen Satz abgeleitet, der sich auf gewisse Systeme von linearen, homogenen Differential- 
gleichungen erster Ordnung bezieht. Ich will einen weit einfacheren Beweis dieses Satzes 
geben und ein paar Bemerkungen beifügen. Es handelt sich um folgendes. 

Es seien beliebig viele, z. B. drei Differentialgleichungen gegeben: 





(de, ER | —( 
di a], T Agdlg 7 Agl, = V, 
dx 
2 
(1) u + bi2, + dan, + bu, = (0, 
dzz 
Fr + 01%, + 69%, + 03 = 0, 


wo die Koeffizienten a,, Qg, Qa, d,,... Funktionen von ti bedeuten. Zu (1) gehört ein 
adjungiertes System 


d 
nn ey =, 
d 

(2) | a bye — ayı =D, 
dy 





ar — 43% — dayg — C3Yy = 0, 


und es hat bekanntlich C. G. J. Jacobi gezeigt !), daß zwischen den Lösungen dieser 
beiden Systeme die Beziehung stattfindet 


(3) %ıyı 4 %oYa + Iays = konst. 
Dies und alles Folgende gilt auch für n Differentialgleichungen, und nur der besseren 
Übersichtlichkeit halber wählen wir n = 3; bei beliebigem n sind alle weiteren Bemer- 
kungen analog zu verallgemeinern. 

Es hat nun Hazzidakis a. a. O. ein neues System dadurch abgeleitet, daß er in (2) 
zu jedem der in der „Diagonale“ stehenden Koeffizienten — a,, — b,, — cz; die Summe 
a, + b, + cz hinzufügte. So ergibt sich 





dz 
> + (ba + 63)2, — bi2; — 43 =0, 
dz 
(4) | z — 422} + (63 + Q,)23 — 0223 = 0), 
dz 
= 43; — ba2, + (aj + b2)2; = 0. 


!) Dieses Journal Bd. 29 (1846), S. 221 oder Ges. Werke, hrsgg. von K. Weierstraß, Bd. 4, Berlin 1886, S. 404. 
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In ziemlich umständlicher Weise zeigte Hazzidakis, daß die Lösungen von (4) 
aus den Unterdeterminanten der durch die partikulären Lösungen x, 23, 23: 2)... .. 
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Zu eir 


Im 90 
einen Satz a 
gleichungen 
geben und « 

Es seiı 


(1) 


wo die Koe 
adjungiertes 


(2) 


und es hat 
beiden Syst 


(3) 
Dies und al 
Übersichtlie 
kungen ana 

Es ha! 
zu jedem d. 
a+b,+c 


(4) 





1) Dieses 
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In ziemlich umständlicher Weise zeigte Hazzidakis, daß die Lösungen von (4) 


!H 


aus den Unterdeterminanten der durch die partikulären Lösungen x, 23, «5: 4. ... 
von (1) gebildeten Determinante 


mn m © 
(5) a’ )- u m u 
„11 ‚I ‚ 


I I I; 


bestehen. Ganz einfach ergibt sich diese Tatsache, indern man zunächst die Gleichungen 
(1) der Reihe nach mit z,, 25, 23 multipliziert und zu den mit x,, X, 2, multiplizierten 
Gleichungen (4) addiert. Man erhält so 


d(2, 2, 7 Lg 2a 2 25) ; 
(6) - Hart) He 


und durch Integration 





Sg Ly2,) = V 


n - 
o , 


(7) ua + usa + a he Jtartbrtene 
wo k eine Konstante bedeutet. 

Um nun die Lösungen der Gleichungen (4) zu finden, trägt man die als bekann! 
vorausgesetzten partikulären Lösungen x1, 23, 23; &/,... von (1) in (7) ein und er- 
hält so 





/ — (ka, +b,-+e,)dt 

Bath at zur ae 
I a |] Er u Z _ (a 1b, -+0,)dt 

(8) a ı+9e 2a+2 z=k e J 1 
#77, ZZ y r — [da +db,-+e,)dt 

\ a 11T ar x 2” Bu k Z e J (a, 2 


wo k’, k’, k'"' willkürliche Konstanten bezeichnen. Bei Auflösung dieser Gleichungen 
nach 2,, 29, 23 hat man z.B. 


/ 


k' 2 T3 
(9\ 21 . (x’ 2" er. — k'' 2 MY E e J (a,+b,+e,)dt 


7 
I SE I 
k 7 ” a 3 


Da nach einer bekannten von Jacobi stammenden Formel !) 


(10) (x . e) =y- u J (a, +b,-+-c,)dt 


wo y eine Konstante bedeutet, hebt sich in (9) der Faktor mit der Basis e beiderseits weg, 
und wenn man y in die Konstanten k’, k’’, k’’’ eingehen läßt, ergeben sich in der Tat die 


ii 


partikulären Lösungen 71, 1,27 von (4) als die drei Determinanten der Matrix 


/ 


(11) | 


Analog erhält man je drei partikuläre Lösungen z, und z, aus den entsprechenden Matrizes. 
Damit ist der Satz von Hazzidakis in einfachster Weise bewiesen. 


IN 


x 0% ©: 


„I 


23 2% X 


Zu den Gleichungen (4) gehört ein adjungiertes System 


dw 
F7 — (bg + C3)wı + Qgw, — Au = 0, 
) dw; | | 
(12) dt bw; — (65 + a,)w + bzw; = 0, 
dw; | | 
ram 27 — (0; + bu); = 0, 





—. 


1) A.a.0. S. 222 bzw. S. 405. 
Journal für Mathematik. Bd. 171. Heft 2. 17 








Dingrldey, Zu emem Satze von Hazzidakis über gewisse Systeme linearer Differentialgleichungen. 


130 


und es finden dann zwischen den Lösungen der vier Systeme von Differentialgleichungen 
(1), (2), (4) und (12) die Beziehungen statt: 


2Yyı 7 Zee ir IYya = c, 
ia ut — | (ay+b,+te,)dt 
Cı 21 ug: 2 192g + Tg2g — ke J ’ 


u) } 
(13) - b) \ - ‘ ! PAOERE ' 
2%] u ZgWs Eee ZgWg =C N) 

k’e* fa +b.+c)dt, 





YıWı — YaWa T YaWz — 


wo c, k, c’, k’ willkürliche Konstanten bezeichnen. 


Sind alle Koeffizienten a,, aa, Qg, d,,... Konstanten, so ist in den vorhergehenden 
Formeln 
(14) eW Jar +b.+e,)dt Bu e* (a, +b, +0) i 
wofür auch e*&+7.+A% gesetzt werden kann, wenn },, Ag, As die Wurzeln der Gleichung 
Cı C5 At, 


bedeuten, die bei Integration der Gleichungen (1) im Falle konstanter Koeffizienten 


bekanntlich eine wichtige Rolle spielt. 
Es sei nochmals daran erinnert, daß sich alle vorstehenden Ergebnisse in sinnge- 
mäßer Weise sofort auf den Fall beliebig vieler Differentialgleichungen ausdehnen lassen. 


Auf besondere Fälle werde hier der Kürze halber nicht eingegangen. 


Darmstadt, den 27. Oktober 1933. 


Eingegangen 28. Oktober 1933. 














